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ABSTRACT 

In the present literature, the storage structures 
whether they are water tanks or silos, are designed as 
per 'elastic' theory with the advent of limit state design 
methods. There is a great need for the limit design of 
such structures. In this thesis cylinderical water 
tanks and storage silos made of reinforced concrete are 
analysed as per the plastic theory based upon the lower 
bound theorem. The pressure distribution considered in 
the case of silos are due to Janssen's and Reimbert's. 

The collapse modes are described and design method is. 
illustrated with an example. 



1. IIMRODUCTIOE- 


1,1 HISTORY OP Civil. EHOINEERIirG 

Civil engineering is a branch, of human activity, 

which has been constantly pursued from times immemorial, 

when people began to adopt the environment for meeting 

their needs. In the ancient times the structures were 

built by artisans. When nomadic tribes inhabited around 

the world, the storage of essentials like water, food etc, 

did not pose any problem, as they never thought of either 

storing or hoarding. With the growth of civilization 

people not only settled at some places but also cultivated 

land to meet their livelihood. This change in lifestyle 

brought into focus the need for storage structures. 

Earthen pots were evolved for purposes of carrying and 

storing water besides storing grains , As the population 

increased, the need for larger storage capacity arose. 

Thus various types of storage structure made of variety 

of materials had come into vogue. The cylindrical water 

(l ) 

tanks and large silos fall under this category , 

The *kno#iow’ of things is known as science and 
‘doing* things is called art. The blend of these two is 
engineering. Civil engineering is the. oldest branch of 
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engineering. Present days civil engineers, structural 
engineers in partictilar, are constantly improving the 
various design procedures, as the behaviour of the materials 
is very well understood with advances in science, 

1.2 IPEiLLIZATION OP STRUCTURES 

An element or an assembly of elements formed into 

a shape to resist a, set of external forces is called a 

( 2 ) 

structure . 

A structure is referred to as a planar or a space 
structure depending on whether or not the centre line of 
the elements forming the struct'ure lie in a single plane. 

All structures are invariably three dimensional, but 
depending on the relative ratios of sides, they are treated 
as one, two or three-dimensional. Whenever it is possible, 
three dimensional structures are reduced to lower dimensional 
categoiy for purposes of analysis and design. Shell, a 
space structure but circular cylindrical axi symmetric ally 
loaded, is idealized as a one dimensional structure for 
analysis and design. 

1.3 STRUOTURAIi ANALYSIS 

Byaluation of the stress resultants, namely, 
internal resisting forces and deformations of a structure 
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caused by external forces* is knoTOi as structural analysis. 

The purpose of structiiral analysis is to study or predict 

the response or behaviour of a loaded structure. The 

entire structural analysis is founded on four basic 

principles, namely the equilibrium of forces, the conditions 

of compatibility, the force-reformation relationship of 

the material used and the knom boundary conditions. If 

all the four conditions are satisfied, the solution is 

Vp) 

exact ^ ,, In this study, the solution found is exact. 

The analysis may be broadly classified as ‘elastic’ 
or ’plastic’, depending on the range of mechanical behaviour 
considered. Elastic analysis predicts the behaviour at 
working loads. The stress at any point in the domain of 
of the structure must not exceed certain permissible working 
stresses. Plastic analysis predicts the behaviour at 
collapse. The stresses at certain points or along the 
lines reach the yield stress of the material. In this 
study, plastic analysis is used to analyse the tanks and 
silos, 

1.4 STRTJOrmAI. PESI&U 

The primary purpose of structural design is force 

' ' t ' ' 

transmission throu^ solid matter with certain design 
objectives such as structural safety, rigidity, serviceability. 
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or geometric, functional and aesthetic requirements. But 
from the engineers 'viewpoint , there are two main criteria 
for the adequacy of a structure which must be satisfied 
in the most economical fashion: Safety, that is sufficient 
strength to resist applied loads without danger of collapse 
and serviceability, that is the ability to carry loads 
without excessive deformation or local distress. 

Ihe strength requirement of the structure may 
be satisfied in several ways. In the 'conventional' 
or 'elastic' design which is also known as 'working stress 
design’, the strength is checked by ensuring that no point 
of the structure is stressed above a value called the 
allowable or working stress i/diich is calculated by using 
factor of safety on the yield stress of the material. The 
actual factor of safety varies from structure to structure 
depending upon the degree of indeterminancy, boundary 
conditions and the nature of loading. 

In the plastic design the total load is investigated 
under vdiich the structure will collapse. The safety of 
the structure is controlled by the choice of a load factor, 
defined as the ratio of the collapse load to the working 
lo^* In this, the behaviour of the structure at working 
load is not taken into acco\mt. 
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In the limit state design, the design load and 
the design stresses are found using partial safety factor 
on load acting on the structure atid the characteristic 
strength of materials. The aim of this design is to 
achieve acceptable probabilities that the structure will 
not become unfit for the use for which it is intended. 

Thus this appears to be the most rational method of 
design to date, 

1.5 MSGHANICAL BEHAVIOUR OF EN&INEERING STRUCTURES 

Structures respond in different ways when subjected 
to external influences depending on their geometric 
characteristics and on the mechanical and other properties 
of the material used. The mechanical properties of a 
material may be represented by determining the functional 
relation between the stresses and strains, 

For most materials encountered in engineering 
practice, it is practically impossible to devise a 
single set of mathematical equations that realistically 
describe the interplay of instantaneous elastic response, 
after-effect and hysteresis, creep and plastic flow. It 
is also a well established fact that unique relations do 
not exist in l^tween stress and strain components in the 
plastic region. The strain depends not only on the 
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final state of stress but also on the loading- history. 

So the stress-strain relations #iich is used in elasticity 
must be replaced by increments of stress and strain in 
plasticity. 

Steel and concrete are the most common materials 
used in structural engineering practice. The strain 
hardening portion of stress-strain curve is neglected 
in the present study. Hence the theoretical computations 
are always on the conservative side of the experimental 
values vhen the idealization is elastic perfectly plastic 
or rigid plastic, 

1.6 IITERATUHE SURVEY 
f 5 ) 

Airy^ derived the equations for pressure 
distribution in shallow bins for granular material based 
on the Coulomb's theory of earth pressure and modified 

(x) 

this for deep bins. Janssen^ also derived equations 
for pressure distribution in deep bins which has been 
followed by Reimbert's deep bin theory^^\ Basically now 
a days the static pressure is calculated by Janssen's 
classical theory or Reimbert's modern deep bin theory. All 
these pressure distribution theory gives static pressure. 

(5 ) 

Brante ,in 1896, made the first measurement on 

( c V 

f ull scale silos. Hoffmann^"^ , in 1916, pointed out that these 
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results showed a lower wall pressure during filling than 
the pressure computed from Janssen's theory, but this 
pressure rises during emptying. 


Cb) (5) 

Subsequent investigations by Bovey'’ and lufft 
in 1904, led to the conclusion that Janssen's formul^a as 
supplemented by Konen in 1896 could, provided that the 
component coefficients were suitably chosen, be considered 
to indicate the wall pressure 'vdien filling and only 
moderate increase were noted during emptying. 


There are two types of flow condition, dynamic 
and nondynamic. The first flow condition gives a large 
pressure increase than the latter flow condition. The 
flow depends upon the filling condition. 

The Reimbert’s equations give a much better 

correlation with experimental data for static and filling 

condition than the traditional Janssen's deep bin theory. 

The emptying condition varies widely with the type of 

( 6 ) 

material and must be considered separately^ 

Save'" and Sawczuk'' have done some vrorh on the 
plastic analysis of tanks and silos for constant hoop 
reinforcement only. Brom the foregoing review it is clear 
that not much work has been done on the plastie analysis 
of tanks /silos with variable hoop reinforcement. In this 
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study tarOcs and silos are analysed for var^ring circumferential 
reinforcement. In the case of silos the pressure distribution 
considered are due to Janssen and Reimbert. 

1.7 SCOPE OP THE THESIS 

The thesis is divided into seven chapters. 

Chapter 1 introduces the fundamental concepts of structural 
analysis and structural design. 

In chapter 2 the analysis, equilibrium equations 
and yield condition of a axis 3 rmmetrically loaded circular 
cylinder shell are presented. 

Chapter 3 deals with the analysis of circular 
cylindrical tanks, while chapter 4 and chapter 5 deal 
with the analysis of silos for Reimbert ’s and Janssen's 
pressure distributions respectively. 

In chapter 6 all the results of chapter 5 have 
been simplified for ultimate load design method and the 
design of a silo is presented. 

In chapter 7 the results obtained in the preceding 
chapters have been critically evaluated and conclusions 
drawn. 
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ESSENTIALS OE AXISYMMETRIClLIiY LOADED OYLINDRIGAL SHELLS 


2.1 DEFINITION 

A sHell is a curved plate. The surface which 

bisects the thickness of plate is called the middle surface. 

If the middle surface and the thickness are specified at 

each point then a shell is completely defined geometrically. 

The ' thickness of a shell is considered to be very small 

compared to the other dimensions of the shell and the 

radii of curvature. Hence a shell can be generalised as 

(9) 

a two dimensional structure' . 

The surface of a shell is obtained by rotating a 
curve of given shape around a fixed axis. For cylindrical 
shell this curve is a straight line. 

2.2 EQUILIBRIUM EQUATIONS 

To establish the equilibrium equation an element 
as shown in Fig, 2i2.1 is considered. The shell of radius R 
is subjected to a radial pressure p and an axial pressure q. 
Due to symmetry the membrane shearing forces N^^q = Nq^ 
vanish in this case and the forces Nq are constant along 
the circumferenee. From sjrmmetry it can be concluded that 
only the forces do not vanish. Also from symmetry it 
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can be concluded that the twisting moments M^q = M0^ 
vanish and the bending moment Mq are constant along the 
circumference. Under such conditions of symmetry three 
(two of rotation and one of translation) of the six 
(three of rotation and three of translation) equations 
of equilibrium of the element are identically satisfied, and 
as such only the remaining three equations should be 
considered. These are obtained by projecting the forces 
on X and z axes and taking the moment about the y axis. 
External forces consist of pressure normal to the surface and 
vertically downward force. These three equations of 

(q) 

equilibrium are^ 



dx R 



dx 


( 2 . 2 . 1 ) 

( 2 . 2 . 2 ) 



(2.2,5) 


Substituting equation (2.2.2) in equation (2.2.1) 
to eliminate from equation (2.2.1) 

+ _ = p (2.2.4) 
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Let Mp and Kp be the plastic moment capacity per unit 
length and axial force capacity per unit length respectively. 
Thus Mp = ’ 



where is yield stress of material and t is the thickness. 
Introducing the non-dimensional parameters 


= 


<1 = 


m^ = 


PR 


a "t 


V 




no = 




R. 


X 


= r 


the equations (2.2.5) and (2.2.4) reduce to 



tR d^x 
4 dx2 



(2.2.5) 


( 2 . 2 . 6 ) 
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2.3 SPaOIPIC POra OP DISSrPATIOlT 

Consider a cylindrical coordinate system x, 9 
and r as shown in Pig, 2.3,1 and also a cylindrical shell 
of length 1 and radius R. Because of the symmetry of 
revolution, circumferential displacements, v, vanish, 
whereas longitudinal displacements, u, and radial 
displacements, w (positive inward) are functions of x 
only. Also due to rotational symmetry, the circumferential 
curvature rate, I'-q, vanishes. There is no circumferential 
displacement. Any point of the shell is displaced in the 
meridian plane in which it is contained. Hence, any two 
neighbouring meridian planes experience no relative 
rotation, and Mq does no - work. So = 0 , The 

radius of the median surface varies from R to R + w, where 
w is the radial displacement of the median surface, the 
circumferential strain, 

/ 

w 

Sq = (-t/2 < z ,< t/2) 

R-z 

must be regarded as constant (and hence K q = 0 ) because 
z (thickness) is negligible with respect to R from the 
very definition of a shell. 

Rates of transversal shear vanishes because it is 
assumed that the material normals remain normal to the 
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deformed medi^ surface. So the dissipation rate is 


h - X % ^6 


(2.3.1) 


where and Nq are generalized stresses, given by 

the bending moment, axial force and circumferential force 
respectively per unit length. 


The corresponding generalized strain rates are 


2* 

dS 

Y = curvature rate = - —r 
^ ax2 


du 


Sx = longitudinal strain rate d= ^ > 

’Jr 

Sq = circumferential strain rate = - - — > (2.3.2) 
y R 


So the equation (2*3*1) will be 


d^w 

D = - % _ + 

dx^ 

using the 'reduced stresses 


du 

dx 


w 
® R 


(2.3.5) 


n^j - Nx/'^p * 

^x ~ > 

ng = Fg/ETp , (2.3.4) 

where the aixial force capacity of the shell element per unit 

length, K^p = a^d the plastic moment capacity of the 

2 

shell element per unit length, Mp = '^y'fc A » vitoere t=thickness 
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of the shell, and <Jy is yield stress of the shell material. 
Therefore, the relation (2.3.3) reduces to 


+ 2 . 2 . 

, , t d w 

B = o„(m^( 

y ^ 4 dx2 


dU 


(2.3.5) 


Hence, the generalised strain rates corresponding to the 
reduced stresses are 



\ 




du 

dx 

♦ 

w 



the relation (2.3,5) is written as 



^ = ‘^y ^ (2.3.6) 

2.4 YIELD OOimiTION 

It is assimed that the concrete is perfectly plastic 
with the yield curve in plane stress as shown in Fig. 2,4.1, 
where and cr are the tensile and compressive yield 
stresses respectively. Here the compressive forces and 
stresses will be regarded as positive, as is usual in 
reinforced concrete practice. Strain rates will be 
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positive when corresponding to contraction. Also 
longitudinal and circumferential directions are principal 
directions. 

It is assumed that the amount of reinforcement 
subjected to compression is small enough to contribute 
in a negligible manner to the compressive strength, and 
that the reinforcement subjected to tension has a yield 
stress Cy and cross-sectional area and A0 per unit of 
length. Because of sjnnmetry of revolution, the strain 
rates are distributed across the shell thickness t as 
shown in Big, 2, 4, 2(a) and 2.4.2(b) for longitudinal and 
circumferential direction respectively, Bor a stationary 
process, the strain rate distribution 'is alike. The 
strain in any layer is 6^ = - z so it is expressed 

in terms of elongation y'V of the middle surface of the 
shell, its curvature ..i<j_and the distance z from the 
middle surface. Due to radial symmetry no circumferential 
curvature occurs. 

-Applying the normality law to the strain rates of 
Big, 2.4.2, it can be concluded that, for the concrete, 
the stress point must be in C (Big. 2.4.I) for 
- ■n t/2 ^ z _< t/2 because e^ < 0 and 62 < 0 and in D 

(Big, 2,4.1) for -t/2 X z _< -T) t/2 because > 0 and 

• ■■ 

82 < 0, During failure the reinforcement yields in both 




^\Csi 






iiii 


m 




I 

' \ 


* 
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axial and circumferential directions, Tb.e generalized 
stresses and Nq can be evaluated for every value 

of the parameter . Elimination of n from the 
expression of and Ng gives the yield surface. The 

reduced generalized stresses are^^^^ 


n 


m 



a' t 

r 


41 

<^2 




The reinforcement ratios are 


tie = Ag/t ; 


and the parameters are 
a 



where t is thickness of the shell, is compressive 

yield stress of concrete, o'j, is tensile yield stress of 
concrete, and yield stress of reinforcement. 

The yield surface is formed of two parabolic 
cylinders given by ; 
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m^(l+a) + + 2n^ + a-1 ) +2p^ \i^ - 4^\x^- 2a = 0 , 

for > 0 (2.4.1) 

p 

-m^(l+a) + 2nx + 2n^ (a_i ) - 2a = 0 , for < 0 (2.4.2) 

This cylinder is bounded by the planes with 
eq.uations 

ng = 1 for compression (2.4.5) 

and nQ = -a-pp-g for tension (2.4,4) 

The tensile strength of concrete is neglected so 
a = 0 . In the case of tank = G but in the case of silos 
Njj is not eq.ual to zero. It is considered in design only. 
Hence the analysis has been restricted to and nQ plane. 
The interaction of the surface with the Hq and m^j for 
given n^ will represent the actual condition of failure. 

The corresponding boundary of the s^e domain in the form 
of a rectangle ABOD is shown in Big. 2 . 4.4 for n^ = 0 . 















3. ANALYSIS OS REINFORCED CONCRETE TANK 
SUBJECTED TO HTOROSTATIC PRESSURE 


5 . 1 GENERAL 

Tanks are generally erected for one of several 
purposes, such as service reservoirs, balancing tanks or 
storage tanks etc. The optimal shape for minimum material 
consimiption is usually a sphere. However when cost of 
erection and others are taken into account, cylindrical 
tanks are the best. In this chapter reinforced concrete 
cylindrical tanks subjected to hydrostatic pressure are 
analysed. The most important aspect besides strength is 
th^, leak proof ness. Thus vdien a fluid storage tank is 
designed, the width of cracks on fluid face must be controlled. 
The porosity of the walls has to be minimised by using a 
rich mix of concrete^^ , 


3,2 EQUILIBRIUM EQUATIONS 


In the case of tank, which is subjected to hydrostatic 
pressure p, the vertical force q is zero. So the dimensionless 


equations (2,3.5) and (2,3.6) will be 


(7) 


dn. 


dx 


= 0 


(5.2v1), 
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— 

4 dx2 " P* (3.2.2i) 

where 

^ t 

P = YX 
Po = Y1 

P"" = t) e = p* ^ 

where 

e = X/D 




24 



+ CG .n0 = 


GO e 
o 


(3.2.4) 


3.3 glassifigatioit Op tauks 

The tanks are classified into three catogeries 
namely shallow, medium, and deep. This olassifleation is 
based upon the mode of failure of the tank at collapse. 
Ibe typical modes of failure are shown in Pig. 3.3.1. ine 
mode Of failure depends upon the value of the non- 
dimensional parameter of the shell CG. 

3.4 method of AHAIiYSIS 


Following are the steps for analysing the tank: 

(i) The distribution of circumferential stress 
resultant Uq is assumed. 

(11) The equation of equilibrium ( 3 . 2 . 4 ) Is solved for 
a given value of the shell parameter 00. 

(lii) The constants of Integration are determined by 

using various boundary and the continuity conditions. 

Thus the distribution of bending moment is determined, 
aare should be taken to see that the shear and moment 
continuous at the critical sections. 


are 
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5.5 SHALLOW TIM 

3.5.1 Stress Distribution 


A shallow tank is also known as a short tank. The 
failure of the tank at collapse is known as mode one 
failure. The distribution of the stress resultants and 
the collapse mechanism are shown in Dig. 3.5,1. The 
distribution of ng is shown in Dig. 3.5.1(c). This is given 
by the equations 


Zone I : ng = anQ for 0 < e < |3 

Zone II: nQ-= anQ + Kn^ (e-p) for P < e < 1 (3.5.1) 

1 - a 

where K = — — 

The boundary and continuity conditions are 


dm^^ 



A I 

®x1 ~ 0 , - 

de 

0 at e = 0 

(3.5.3) 

and 

dm^l 

m^j “I ~ m^ 2 Sind — — 

de 

dnn. p 

= — at e = P 
de 

(3.5.4) 


3,5, 2 0 < e < P (Zone l) 


In this zone, which is from 0 to p, the circumf erential 
reinforcement is constant and equal to ocnQ. So the equilibrium 
equation ( 3 . 2 . 4 ) reduces to 



CO (p*e - an^^) 


(3.5.5) 






iN^S®!i;Ki;PSS 

.:y‘;.,- ■’. 

,'... •■*■ ■■-•.’V.-v 
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Integrating once 
dm. 


de 


-n* 2 

•V "I P ^ ® 

= OG (-2 — 


~ + A 


1 


(3.5.6) 


Integrating once again 

.3 ^ 

/-o 

^“7 


Po e5 an^e 


-|— ) + A^e + 


(3.5.7) 


wiiere and are the constants of integration. 
Prom bomidary condition (3,5,3), one gets 


Ai = B^ 


= 0 


(3.5.8) 


So equations (3.5.6) and (3.5.7) will now be 


dm^ji 

, ^"0 e 

- an^e) 


de 

= 00 (-|— 

( 3 . 5 . 9 ) 


,Po®^ 

an e^ 


“x1 

= oc ( ° - 


(3.5.10) 


3.5.3 P < e < 1 (Zone II) 


In this zone, which is from p to 1, the circumferential 
reinforcement is -varying linearly from an^ to n^. So the 
equilibrium equation (3.2.4) will reduce to 

d^ 

2 ' ^ = ^ ® - o^i^o ’■ ^^0 ( e - P )) (3.5.11) 


First integration gives 
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dm^2 , Poe^ 

= CG ( - an^e - Kn^ ( ~ - pe ))+ (3.5.12) 


Second integration gives 



GO 



tr / V 

~ (3.5.13) 


where and Bg are constants of integration. 


Prom the continuity conditions (3.5.4) 

r... o2 


Ao = - CG 


Kn^P' 


Kn^p 


(3.5.14(a)) 

(3.5.14(b)) 


So equations (3.5.12) and (3.5.13) will be 

Kiln , ,2, 

~ - CO (-1 an^e - (e-P) ) 


- '‘“■0® - 


CG (- 


an^e' 


(e-P)b 


Ihere is one more boundary condition 


“^1 = ®o e = I 

Using this in equation (3.5.16) 

GO 

“o = — (p* - 30tHn - Kiig 

and after simplification one gets 


(3.5.15) 

(3.5.16) 

(3.5.17) 

(3.5.18) 
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6nio ■ 2 ; 


(3.5.19) 


.By using the equations (3.5.10) and (3.5.16) , 
the bending moment distribution will be as shown in 
Jig. 3.5.1(b). 


3 . 5.4 Gritical Condition 

As the Ii/R ratio of the shell increases, the value 

of GO will also increase and from equation (3.5.19) it is 

clear that the value of p* will decrease for increase in 

the value of GO , From equation (3.5.16) one can infer that 

the increase in 00 will cause decrease in m^ and for some 

value of e , and GO , m^j will reach the maximum negative 

value^^^\ So for certain value of e = e^ 

dm,. 

m^ = - m^ and =0 (3.5.20) 

when ©Q i< P 

dm^ji 

= - io ana = 0 at e = 60 


This gives 

SuzIq 


®o = 


and GGR = 


6 m, 


3anoe^ - p*e^ 


(3.5.21) 

(3.5.22) 
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where COR is the critical -^alue of GO. 


When e > P 

0 - 


m ^2 = - Dirv and 


o 


= 0 at 


de 


e = e. 


This gives 

2 9 1/2 

(anQ-Kn^P) + ((anQ-En^P) + Kn^Cp^ - KnQ)p^) 


®o = 


OOR = 


(v* - aig) 

6 mrv 


(3.5.23) 


(e^(KnQ-p->) + 3noeo(a>Kp) + 3EnoP ^eQ-Kn^^p 


3.5.5 Special Case 


let = n^ = 1 


Then p* = 6/CG + 3 a 


Ror zone 1 , 0 < v, _< p 


(3.5.24) 


®o = 2 q;/p* 


OOR = 6/(5ae^ - K 


(3.5.25) 


Ror zone 2, P .< e ^ 1 


(a-KP) + ((a_Kp)‘=^ + K(p* - K)p ) 


2 . 1/2 


(Pj - E) 
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OCR = 


(e^(K-p*) + 56^ (a-KP) + - KP^) 


(3.5.26) 


Let a = P » 1 


Then K = 0 


p* = 3 ( 1 + 2/GC) 


®o = 2/p; 


OOR = 




(3.5.27) 


3.5.6 Results 


The distribution of is dependent upon tlie values 
of a and p. Ror a given set of a and p values the critical 
shell parameter value COR, the location of hinge e^^, and 
the collapse pressure p* are shown in Table 3.1 for mQ=mQ=nQ 


TABLE 5.1 : OCR ROR RIRST MOLE OF RAILURE CR TANK AT 


®o = “o = “lo = 


a 

P 

COR 

®o 


0.25 

0.25 

118.5 

0.52999 

1.22251 

0.50 

0.25 

39.6 

0.58205 

1.93500 

0.50 

0.50 

38.0 

0.56509 

1 .78289 

0.75 

0.50 

25.3 

0.58631 

2.57057 

1.00 

1.00 

17.1 

0.59685 

3.35088 
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3.6 MEDIUM TAM 

3.6.1 Stress Distribution 


A medium tank is suck that there is hoop compression 
at the free end. This type of failure is designated as 
mode two faiilure and this is termed as ’total failure’ or 
total collapse. The distribution of stress resultants 
and collapse mechanism are shown in Dig. 3.6.1. The 
distribution of n^ is as shown in Dig. 3.6.1(b). There 
are three zones: 


zone 

I 

* 

Zone 

II 

: n, 

Zone 

III 

: n, 

The 

boundary i 


= . 

0 » 

“x1 

= 


ro 

= 

“x3’ 


g = -'^nQ I or ^ \ e^ 

g = ocn^ for e-| ^ e ^ P 
g = an^ + K (e-P) n^^ for P ^ e ^ 1 


^x1 

de 

^x1 

de 


= 0 at e = 0 
dm„p 

= — at e = e. 

de 

= ^ at e = P 
de 


( 3 . 6 . 1 ) 

( 3 . 6 . 2 ) 

(3.6.3) 

(3.6.4) 


The continuity conditions at e = e 2 depend whether or not 

ep > p 
< 

If ep ^ P 
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1^2 = - 


dm^2 


~ 0 at 


© — 0 / 


and if 


e, > P 


“xS ”o > 


an^3 


— 0 at 


e = e. 


(5.6.5) 


m at e = 1 


( 3 . 6 . 6 ) 


3,6,2 (Zone l) 


In this zone, which is from 0 to e^ , the circumferential' 
stress IS negative and reinforcement provided is constant 

and equal to an^ , So the equilibrium equation (5.2,4) 
reduces to 




- GO (p*e + aiiQ) 


Integrating once 


(3.6.7) 


OG ( + an^e) + 

I3i'feeg2?a"fcing Qnce again 


“x1 = CO (- 


arioe- 


-A-j© + B-j 


where and are constants of integration, 
i'rom boundary conditions (3.6.2) 


A. 3 j 


(3.6.8) 


(3.6.9) 


0 


( 3 . 6 . 10 ) 
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So eq.uations (3.6,8) and (3.6.9) will be 
__£i — ^ QQ + an^e) 

de 

p*e^ an^e 

and iiLjj.-,= GO ( ) 


( 3 . 6 . 11 ) 


( 3 . 6 . 12 ) 


3.6.3 < e < P (Zone II) 

In "blais zone^ wiiich is from e^ top ,th.e circumferential 
stress is positive and reinforcement provided is constant 
equal to an^ from e^ to p. So the equilibrium equation( 3. 2, 4 ) 
will be 


V2 2 N 

= C (p*e - an^) 


(3.6.13) 


From first integration t one gets 


JP^e 

= GO (-| an^e) + A 2 


(3.6.14) 


From second integration , one gets 


:2 = GO ( 


■) + A^e + B2 


where A 2 and B 2 are constants of integration. 
Using the continuity conditions (3.6,3) 


±2 = GO ae.j (n^ + eLq) 


(5.6.15) 


GGae. 


®2 




2 


( 3 . 6 . 16 ) 
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So equations ( 3 . 6 . I4) and (3.6.I5) reduce to 


dm. 


•x 2 


de 

and 


P*e 

- CJG ( an^e + ae^ (n^ + n^)) (3.6.17) 


ocii^e 


%2 CO (-^ . -j£_+ ae^eCn,, + H„) - 


( 3 . 6 . 18 ) 


5 . 6,4 P < e _< 1 (Zone III) 


this zone, ■which is from P to 1 , the circumferential 
stress is positive and reinforcement provided is linearly 
varying from an^ to n^ from P to 1 , So the equilibrium 
equation (3.2,4) reduces to 


,2 

d m 


de 


■x 3 

2 - an^ - Kn^ (e-p)) 


(3.6.19) 


J'irs't in’tegra'tion ogives 

^%3 _ . e2 

— - CG (_ - an^e - Kn^( ~ - Pe))+ (3.6.20) 


Sscond. iirfc 0 § 3 ?a,"bio 2 i gives 


,P^e- an.e‘ 


R 2 

e pe 


m_^ = GO (-S .. ,Zo : ,, . f 

2 ^”■ 0 ^ g ” 2 ))+ A^e + 


where and are constants of integration. 

Prom the continuity conditions (3.6.4) 

A3 = 00 (ae/„„ + 5„) ) 


( 3 . 6 . 21 ) 
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Kn.p^ 

B = CG ( — 

5 6 

So equations (3. 6 . 20) and (3.6.21) reduce to 


ae-i 

-7- 


(3.6.22) 


j%3 

de 

^ 2 

o^n^e - Kn^ 

“x3 

p*e^ 

= 00 (-2^ 

^2 

anoS 

2 


Kn^ep^ 

2 

-h ^ - 

6 


KnoP' 


Pe2, 


ae 


(3.6.23) 

+ ae-,e (n^ + 5^) 

+ ^o^) (3.6.24) 


Using boundary condition (3.6*6) in equation ( 3 . 6 . 24 ) 
one gets 

6mo 

J'S = ST ^ + KHo (1-P)5 . (5^ + 

( 3 . 6 . 25 ) 

Prom the first condition of ( 3 . 6 . 5 ) 

^ _2(an„e2 - ae^ (n„ + S^)) 

2 

' 2 ^ 




(e?) 


(3.6.26) 


and 


PS 


_ (-6mQ/tlC+ 5anoe2 - Sae^e^ (n^ + n^) + 5 ae^(nQ 


n. 


,)) 


(e? ) 


(3.6.27) 


Prom the second condition of ( 3 . 6 . 5 ) 

p. = i!“o®2 - 2“ei(no + 5^) + Ka„(e 2 -P)^) , 

^0 - -• - ^ (3.6.28) 

( ef ) 


and 
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(-65^/00 + 5anQe| - 6ae^e2(nQ+no)+2ae^(rLjj+nQ)+KnQ(e2'-P)^) 

Po=“~^ — — 

(5.6.29) 

By using equations (3.6.12), (3.6.15) and ( 3 . 6 . 24 ), the 
bending moment distribution will be as shown in Pi'g. 3.6,1 (a). 


3.6,5 Critical Condition 


As the value of CC increases the maximum value of 
positive moment increases at certain value of CG at e = e^. 
This value of CG is critical. This is clear from the 
equation (3.6.12), The value of CCR will also depend on 
whether ©3 P . The value of e-j is confined upto O.IL 
of the upper portion, so there is no point in considering 

P < ©I* 

So at e = e^ 

m^ = + m^ and 0 ( 3 . 6 . 30 ) 


when 


e, i P 


- e. 


o 


(3.6.31) 


CCR = 


6mr 


(pS®5 - 3anQe^ + eae^e^Cn^+no) - 3ae^(nQ+nQ)) 

( 3 . 6 . 32 ) 



40 


and when 2 P 


2ng( g - KP ) 
(p* - Kno) 


- e. 


(3.6.33) 


OCR = 


6mr 


(p*e|-3anQe|+6ae^e^(nQ+nQ)-3ae^(nQ+nQ)-KnQ(e^-B)^) 

(3.6.34) 


3 . 6.6 Special Case 


When m = 5 q = n = Sq = 1 


p* = 6/CO + 3a + K (1-p)^ - 6ae^(2-e^) (3.6.35) 

when 62 2 P then equations (3.6.26)and ( 3 . 6 . 27 ) will be 

-6/00 + 3ae| - 12ae^e^ + 6aef (3.6.36) 

o ^3 


®2 - 


®2 

2 1/2 
a + ( a ~ 4aPo®i ) ' 

n 


(3.6.37 ) 


When 62 2 P then equations (3^6.28) and (3.6.29) reduce to 


p* 


0. 


^/OG 12 o 6 e.j 62 6(xe^ g 38) 


- ± .( (a-KP)-(Po-^)(4ae^- ^ 


1/2 


(p; - E) 


When < P equations (3.6.31) and ( 3 . 6 . 32 ) will be 



CGR 


6 


(3.6. 4-1 ) 


- 3ae^ + I2ae^e^ - 6ae^) 

When 2 P then equations (3.6,33) and (3.6.34) will be 


®3 = 

2( a - Kp ) 

®2 

PS -K 

(3.6.42) 

GGR = 

6 

(pSe|-3ae^+1 2ae^e^-6ae^-K(e5-P )^) 

(5.6.43) 

When 

II 

CCL 

(I 


P* 

■*^0 

= 6/CG + 6e^ - 12e^ + 3 

(3.6.44) 

®2 

2 ± r 1- 4PSei 

P* 

(3.6.45) 


— 0 



2 

= - ©o 

p* 2 

0 

(3.6.46) 

OCR = 

6 

(3.6.47) 


(P*e| - 3e^ + 12e^e^ - 6e^) 


3.6»7 ''Qritical Values 

Ihe distribution of n^ is dependent upon the values 
of a and p, For a given set of a and P -^ues.the critical 
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shell parameter -value COR, location of hinges 62 and e^, 
the range of hoop compression from 0 to e^ , and the collapse 
pressure p* are shown in Table 3.2 for m^ = m^ = n^ = n^ = 1 . 

3.2 : GCR ROR SECOm) MODE OE EAILURE OP IMK EOR 

“'o = "^o = ^ 


a 


CGR 

®1 ■ 

®2 

®3 


0.25 

0.25 


— 

— 

... 


0.50 

0.25 

250.15 

0.075 

0.736 

0.2071 

1.37369 

0.50 

0.50 

- 

- 

- 

— 

— 

0.75 

0.50 

147.20 

0.0862 

0.7377 

0.16436 

I.6IO233 

1.00 

1.00 

93.10 

0.0890 

0.7446 

0.2349 

,2,043000 


As the value of CG approaches CGR the medium shell 
becomes deep shell and failure mode changes from second 
mode to third mode, 

5.7 D1E2 TAIK 

3.7.1 Stress Distribution 

This tank is also knom as long tank. The failure 
of the tank at collapse is knovoi as mode three failure. The 
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distribution of the stress resultants and the collapse 
mechanism are shown in I^ig, 3.7.1. This mode of collapse 
is only partial. The upper portion of the tank remains 
intaot. The distribution of Uq is shown in Eig, 3.7.1(b). 
This is given by the equations 


Zone I : ^0 = ® ^ P 

Zone II : Ug = an^ + K (e-P) for P ^ e ^ 1 (3.7.1) 


The boundary and the continuity conditions are 
dnu-i 

“ " ■ fore^<P 


m 


x1 


m, 


o ’ 


%2 = “o » 

= “ 3:2 » 


de 

dmx2 


= 0 at 


6 = 6; 


de 


- 0 at 6 = 6^ for 2 P (3.7.2) 
at 6 = p for P 


de 

dm^i _ dm^2 


de 


(3.7.3) 

The continuity conditions at e = e^ depends whether 

e. > p 

5 ^ 


Eor eg ^ P 


_ 


= 0 at e = 


er 


«S:1 = - “o " ae 

and for 

^ >5 >>■ - ' 

*^2 

“x 2 ~ ^ ®o * ~ = 0 at e » Sc 

de ' 

and 


( 3 . 7 . 4 ) 
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™x2 ~ e = 1 


( 3 . 7 . 5 ) 


3.7.2 < e X p (Zone l) ' 

In tills zone, circiffiif erential stress is positive 
and reinforcement is constant throughout, from e^ to p, 
is an , So equation (3.2.4) vdll be 


■^1 


de 


^— = oc ( p*e - an^) 


Integrating once 


ae 


( Pj ^ - a n^e ) + 


( 3 . 7 . 6 ) 


( 3 . 7 . 7 ) 


anoe 


) + A.e + B. 

4 4 


Integrating once again 

, pS®^ 

^1 = 00 ( 2 

A^ and are constants of integration 

Prom the continuity conditions (3.7.2) 
p*e^ 

A4 = - 00 ( - an^e^) 

p*e? an_e^ ^ 

= m, + 00 ( -|-i f—) 


(3.7.8) 


(5.7.9) 


Substituting these values in equations (3.7.7) and (3.7,8) 


one gets, 
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dm. 


de 


x1 _ 


= OG (■ 


2 2 
p*e P*e7 

— - - an^e - + an^e^) 


(3.7*10) 


^ an^e^ P*e?e p*e? 

= GO - —2 -2-i + an e.e + -2_±- _ .--0-i)+ 

6 2 2 0 4-^ 2 o 


3 2 
(3-7.11) 


3.7.3 P < e < 1 (Zone II) 


In this zone circumferential stress is positive 
from p to 1 , and the reinforcement provided, in this 
region is linearly varying from an^ to n^. So the equation 
( 3 . 7 . 4 ) vd.ll reduce to 

d^ 2 

■ — “2 — ^ ^^o® " “’^o " ^o (3.7.12) 

de 

Pirst integration give 
'^‘^2 


de 


= <3C ( -2— - an^e - Kn„( — - Pe)) + A, 


(3.7.13) 


Second integration gives 
x2 = « (-j 2— -Eno (• 


m, 


Pe^ 

- - — )) + k^e + B 5 


( 3 . 7 . 14 ) 


vshere and are constants of integration. 


Prom the continuity condition ( 5 . 7 . 5 ) 
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Pjef Ka 

Aj = CC (cxn^e^ - 2 ) 


B 5 = 


+ OG ( 




Pq Q 4 _ °^^o^4 + ) 


( 3 . 7 . 15 ) 

(3.7.16) 


Substituting these values in equations (3.7.13) and 
( 3 . 7 . 14 ) one gets 


dm 


•x2 , , 2 2. 


de 


= GG ( _(e2 2) ^ an^(e-e,) - (e-P)'") 


(3.7.17) 


and 

m^2 = 00 (^ ( 


( 1! - e^e !i- ) . 1°. (e e )2 

2^3 ®* 4 ® ^ 3 ^ 2 ^ V 


Ku 


a(e-P)5) + 


m. 


(3.7.18) 


Using boundary condition (3.7.5). 


p* = 


3anj^(l-e^)^ + En„(l-P)^ 
(l+2e^) (l-e^)‘ 


(3.7.19) 


3.7.4 < P (Case I ) 

For p, using the continuity conditions (3.7.4) 
in equation ( 3 . 7 . 11 ), one gets 

* ■ 

Pq = 2anQ/(e^+e5) (3.7.20) 

and 
OG = 


e(m^ + Sq) 


{ccn^{3e^-6e^e^+3ep - pj (e^-3e^e^ + 2e^ ) ) 


(3.7.21) 
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when P ^ > using tlie continuity conditions (3-7.4) 

in eq^uation (3. 7. 18), one gets 

,2 

(3.7.22) 


2 anQ(e^-ep + Kno(e^-P)' 

Ir 


( 2 2^ 
(e^ - e^) 


6(mo+mo) 


OC = 


(3.7.23) 


n 2 . 


( 3anQ ( e^-e ^) ^+KnQ ( e^-p ) ( e|-3e^e^+2e|) ) 

3.7,5 P i (Case II ) 


lor this condition equilibrium equation (3.2.4) is 


dV 


■x3 _ 


de' 


= CO (p* e - an^ - 


(3.7.24) 


First integration gives 


<im 


x3 


. e^ 


de 


= CG(-~ - an^e - Kno( ^ - Pe)) + Ag (3.7.25) 


Second integration gives 

P*e^ oLne^ 

■^3 = 0° c V S- - V - 


(3.7.26) 


where Ig ahd Bg are constants, of integration. 
Using the continuity conditions (3.7.2), one gets 
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“^6 “ 


/ ^o®4 
CC ( 




Kn^ (-^ - Pe^)) (3-7.27) 


and 


= m.— CO 
6 o 


P*0«I 

'■ - 3 2 


e? Pe| 
Knot ^ - -J- 


(3.7.28) 


Brom the boundary conditions (3.7.5) 


- ( l-e^j^+Kn^C 1-.3P+6Pe^-3pe^~3e^4-2ep 


pj 


(1- 5e|+2e|) 


and from the continuity condition (3.7.4) 


(3.7.29) 


®5 ®4 = 2(a-KP)/(p*-K) (3.7.30) 

6(mo+5o^ 

Oc = — 

( 3ano ( e^ -e ^ ) ^+En^ ( e|-3P e^-3e^e^+6§ e^e^+2el-3§ ef ) 

-PS(e|-3e^e^ + 2e|)) (3.7.31) 

3.7.6 Special Case 


Let = 5^ =1 

For ^ p equation (3.5.19) will be 

pj = (E(l-p)^+3a(l-e^)^)/((l+2e^)(l-e^)^) (3.7.32) 

When < p, equation (3.7.20) and (3.7.21) will be 
P* = 2a/(e5+e^) (3.7.33) 
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and 

00 = 12/(3a(e^-e^)%*Ce|-5e|e5+2e^)) (3.7.34) 

when p _< , equations (3.7.22) and i3*l»23) reduce to 

P* = (2a(e5-e^)+K(e5-P)^)/(e5-e^) (3.7.35) 

and 

00 = 12/(3a(e5-e^)^+K(e5-e^)5-pJ(e|-3e^e5+2e|)) (3.7.36) 

when 2 equations (3,7.29) , (3.7.30) and 

(3.7.31) will be 

(l+2e^) (l-e^)^ 

12 

GO = — - — 

(3a(e^-e^)^+K(e|-3Pe|-3e|e^+6pe^e^ + 2e^ 


-3§e^) - p*(e|-3e^e^+2e^)) (3.7.58) 

and e^ + e^ = (2(a_E;p))/(p*-K) (3.7.39) 

idien a = p = 1 

P; =3 / (l+2e^) (3.7.40) 

00 = 3(l+2e^) / (1-ep^ (3.7,41) 
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(5-7.42) 


By using tlie equations (3.7.11) and (3.7.18) or 
(3.7.26) the bending moment distribution will be as shown 
in Big. 3.7.1(b). 

3,8 COITCLUSION 

The graph is drawn for the collapse pressure p* 

vs shell parameter CC as shown in Big. 3.8.1 for various 
values of a and P; Brom this graph it is clear that 62 aud 
e^ are always greater than p and e^ and e^ are always less 
than p, for P less than or equal to 0.5 and a upto 0 , 75 . 
Hence the curves are smooth and only changes in slope 
are noticed at critical points. Brom this figure it is 
also clear that the negative circumferential stress is only 
confined upto 10 percent length of the top portion of the 
tank. Bor any given tank the collapse pressure and 
corresponding mode of failure can be determined by 
calcialating the shell parameter GC, 











liiPI 


different oc 





4. ANALYSIS CfP REIIC’ORCED CONCRETE SILOS BASED ON 
THE REIMBERT'S PRESSURE DISTRIBUTION THEORY 

4.1 G-INERAL 

Large size containers to store grain, cement, coal 

etc. are temed as ’bins*. A container ■vdiich. is sHallow 

is known as a 'btmker*. In this case the plane of rupture 

meets the grain surface at top. So hy proportions, a bin 

would be called a bunker if its height L is less than 1,5 

times the diameter . D for circular shapes and if height L 

is larger than 1.5 times the diameter D it is called a silo 

In the case of silos the plane of rupture meets at the 

C 1 1 ) 

opposite side wall , 

The shape and size of the silos are determined by 
the engineer by taking the site, storage or any other 
service condition. In the silos the width of crack 
allowable is more than that of the tanks, A niimber 

(12) 

of pressure distributions are given in literature . 

In this thesis ,jcy 11 nderical. silos subjected to Reimbert’s 
and Janssen’s pressure distribution are analysed. In 
this chapter plastic analysis of cylindrical silos 
subjected to Reimbert’s pressure distribution is presented. 
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4.2 PRESSURE DISTRIBUTION 

The material filled in cylindrical silos, in ^ich 
emptying hole situated at the centre of the bottom, is 
held in equilibrium by the reaction from the sloping bottom 
and the frictional force on the vertical wall, as shown 
in Pig, 4.2,1^^^^. 

The lateral unit pressure at depth X from top is 
P = ( ■> - ( 1 + X/G)‘^) (4.2.1) 

Por circular silos 

W = YI>/(4H') (4.2.2) 

md 0 = D / (4n'l!:) - lig/3 (4.2.3) 

where 

p* = tan 0' , coefficient of friction between stored material 
and silo walls, and 
k = (l-sin 0)/ (1+sin 0) 

0t = an^e of friction between material stored and silo wall, 
hg == de^th of the heap, 

0 = angle of internal friction of stored material, 

D = dismeter of silo , ^ ^ ^ ^ ^ ' 

R^ = hydraulic radius , and 
Y = unit weii^t of stored material. 
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For deep silos, li„ is taken as zero. Substituting 

o 

these values in equations (4,2,2) and (4.2.3) » one gets, 


Pmax =YRi/1^’ 
0 = 


So equation (4, 2 . 1 ) reduces to 
P = 


YR^ XP'lc 

—4. (i_(i ^ — 


Rh 


-2 

) ) 


(4.2.4) 

( 4 . 2 . 5 ) 


(4.2,6) 


The vertical static unit pressure at depth X below the 
surface of the stored material is 

q = Y (X ( X/G + 1)“I hg/ 3 ) ( 4 . 2 . 7 ) 


Substituting the value-, of 0 and taking hg as zero, 
equation ( 4 , 2 , 7 ) reduces to 

q = Y (X (Xp’k/Tl^ + 1 )“^) (4.2.8) 

The vertical frictional force on a unit width of wall at 
depth X is 

y = ( YX - q ) R^ (4.2.9) 

Static unit pressure normal to a surface inclined 
at an angle 0 to the horizontal at depth X below surface 

(12) 

of stored material is 

q^ = P sin^Q + q cos^B > (4.2,10) 
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4. 3 EQUILIBRIUM EQUATIONS 


Ihe eq^uilibrium eq^uations for a silo in non- 
dimensional form are given by the equations (2.3.5) and 
(2,5.6), Unlike water tanks there is a vertical pressure 
along the walls besides the normal horizontal. The 
^quilibrim equation for the normal pressure is given by 


■feR d^ 


+ ng 



(4.3.1) 



(4.3.2) 

(4.3.3) 

(4.3.4) 

(4.3.5) 


e = X/E, (4.3.6) 

where 

Y = unit wei^t of stored material, 

R^ = hydraulic radius , 

U* = coefficient of friction between stored material 
silo walls, 

T = nondimensional friction factor, 
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r 

t 

OG 

X 

I 


crushing stress of concrete in compression, 

= thickness of the shell wall, 

= non-dimensional geometrical parameter of shells, 
= depth from top and 

= total depth of silos. 


Substituting equations (4.5.2) and (4.3.5) in equation 
(4 .3.1), one gets 


2' + CJO ne = GO p* ( i_ (1 + el)“^) 

Q.© 


(4.3.7) 


4.4 SHORT sno 

4,4.1 Stress Distribution 

A shallow silo is also known as a short silo. The 
failure of the silo at collapse is known as mode one 
f ailure as in tank. The distribution of the stress 
resultants and the collapse mechanism are shown in Fig, 4,4.1 
The distribution of n^ is as shown in Fig. 4 . 4.1 (c ) , This 
is given by the equations 

Zone I : n^ = an^ for 0 ^ e _< p 

Zone II I Uq = an^ + Kn^j (e-p) for P X e < 1 

Tflhere 


K ' = 


1-a 
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The boundary atid the continuity conditions are 


dm. 




“x 1 - ”22 > 


x 1 


de 

dm 


x 1 


- 0 at e = 0 

dnix2 


and m^^ = m^ at e = 1 


4,4,2 0 _< e ^ P (Zone l) 


de 


at e = p 


(4.4.2) 

(4.4.3) 

(4.4.4) 


In this zone, which is from 0 to p , the circumferential 
reinforcement is constant and equal to an^. So the 
equilibrium equation (4.3.5) will be 


d%ji 


-2 


— = GG ( p; ( 1,.(1 + eT) _ an ) 

de'^ ° o' 

Integrating once 

, -,-1 , . . ..-1 


(4. 4 . 5 ) 


de 


= GG (p* (e + T" ( 1 + eT)“ ) - an^e) + 

(4.4.6) 


Integrating once again 


• ^ (P^ ( ^ T ln(l+eT)"')- )+A^e+B| 

(4.4.7) 

where and are constants of inte,gration, 
i'rom the boundary conditions (4,4.2) 


A^ = - GO p* T 


-1 


(4.4.8) 
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= 0 


(4. 4-. 9) 


So equations ( 4 . 4 » 6 ) and (4. 4 * 7 ) reduce to 


dm 


•x1 


- CG (p*(e+T ^( 1 +eT) T’’”' ) - anQe) (4,4,10) 


de 


and 


o 4 

(P*^ ® + ^ ln(l+e!I!) - T“'e) - — ) 


(4.4.11) 


4.4.3 P ^ 1 (Zone II ) 


In this zone, which is from p to 1, the circumferential 
reinforcement is varying linearly from an^ to n^. So the 
equilibrium equation (4.3.5) becomes 


de^ 


- CC (p*(l - (1 + el)"^) - an^ - Kn (e-p)) 


(4.4.12) 


I* irst integration gives 


dm. 


X2 


de 


CO (p* (e + T ”*( 1 + 61 )”^)- anQe-KnQ( Pe))+ A 2 


Second integration gives 


(4.4.13) 


2 '=5 ? 

/ /® -2_ pe^ 

= ac(p*(^ I ^lad+eT))- -f— - i- )) 


^2® B 2 


6 2 

(4.4.14) 


wher© A2 ^id B2 3re constants of integration. 
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From the continuity conditions (4.4,3), one gets 

. P^n^ 

A2 = - CG (V* T"' + — ) (4.4.15) 

GC KUq 

= ( 4 . 4 . 16 ) 


Substituting the values of A2 and B2 in equations (4,4.13) 
and (4.4,14) , one gets 


y2 -I 11 2 

~ CO (p*(e+T”' (l+eT)~ - 1 “ )-an e - —r- (e-P) ) 

de ° 0 2 

( 4 . 4 . 17 ) 

and 2 2 

o 1 o^n„e Kri ■3; 

m^2 = (P*(— + T‘''^n(l+e!P)-T"' e) - -^(e-P)^) 

( 4 . 4 . 18 ) 

From boundary condition (4.4,4) 

6mo/GD + 3ano+KnQ(l-p)5 

p* = — 2 ^ — 2 2 .. (4.4.19) 

3+6T”^ln(l+T) - ST"’' 


By using the equations (4.4.II) and ( 4 . 4 . 18 ) , the bending 
moment distribution will be as shown in Fig. 4.4.1 (b). 

4.4.4 Critical Condition 


As the ratio of the shell increases, the values of 
GC as well as 1 increase . If L/k ratio remains same, but 
both 1 and R increase the value of OG increases 2 ' 
remains constant. As CG increases decreases, sis is evident 



63 


from eq.uation (4.4,19). As p* decreases, m^^ and 
also decrease which is evident from equations (4.4,11) and 
(4.4«18). But for certain value of OC = COR and a particular 
e = Sq, or m^2 'wiH reach the maximum negative value 

dm-1 

HLy = - 5- and = 0 (4. 4. 20) 

^ ^ de 



fhen 
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P* = 


6/GG + 5a + K(l_p)^ 


° 5(l+2T"^ln (1+T)-2T“‘* ) 


(4.4.25) 


when, < P 

A — ‘ 


and 


P* = 


GCR= 


a (1+6^1) 


a e^-p* (eQ+ZE-^lnd+egT ) -2T-'' ) 


(4.4.26) 


(4.4.27) 


when p _< e _< 1 


= 


2ae^ + K (e^ - p)^ 
2(e^ +!E"'^(l+e^!r)””‘- (D’*”' 


and 


(4.4.28) 


GGR = 


( 3ae^+K ( e^-P ) ^-3 p^ ( eo+2T'“^ln( 1 +eoT )-2!I}"'' e^ ) 

(4.4.29) 

^en a = p = 1 

2 ■ 


p* = . , ,,. 

0G(l+2T"'^ln(l+!r)-2T” ;) 


(4.4.50) 


= 


5(p* - 1) 


(4.4.31) 


GGR = —y " ^ ' ' ' - ^ — ■ 7 " ' ' i ^ ■ 

eo-P*(eo+2!r”nn(1+epT )..2T"^ e^) 


(4.4.32) 
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4.4.6 Results 

ITie distribution of Uq is dependent upon the values 
of a and P. S' or a given set of a and P values and for 
prescribed T, the critical shell parameter value OCR, the 
location of hinge e^, and the collapse pressure p* are 
shown, in Table 4."1» 4.2 and 4.3 for mQ = m = n^ = 1 and 
for T equal to 1 , 2,5 and 4 respectively, 

TABLE 4.1 : OCR E OR FIRST MODE OP FAIIBRB OP SILOS AT T=1 
. AML = 5o = Uq = 1 


a 

p 

OCR 

®o 


0. 25 

0.25 

267.5 

• 0.338 

1,0300 

0.50 

0.25 

66.7 

0.506 

1.6140 

0.50 

0.50 

58.5 

0.504 

1.4897 

0.75 

0.50 

35.0 

0.539 

2.1463 

• 

o 

o 

1.00 

24.6 

0.549 

2.7490 
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TABLE 4.2 : GCR EOR EIRST MODE OP EAILURE OE SILOS AT 
T = 2.5 AND hIq = iiIq = = 1 


a 

P 

OCR 


pr 

0.25 

0.25 

515.0 

0.249 

0.656 

0.50 

0.25 

111.8 

0.424 

1 .018 

0.50 

0.50 

88.5 

0,455 

0.939 

0.75 

0.50 

50.2 

0.500 

1.350 

1.00 

1.00 

34.8 

0.518 

1.760 


TABLE 4.3 

: OCR FOR FIRST 

T=4 AMD mQ=5Q= 

MODE CP PAIHJRE CP SILOS AI 

="0=-' 

a 

P 

OCR 

®0 


0.25 

0.25 

819.0 

0,201 

0.560 

0.50 

0.25 

163.0 

0,362 

0,^4 

0.50 

0.50 

119.0 

0.421 

0,796 

0.75 

0.50 

65.0 

0.477 

1.144 

1.00 

1.00 

44.2 

0.494 

1 .490 
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I’rom the forgoing Tables 4-1 to 4.3 , it is 
clear that as the value of T increases the range of OCR 
for first mohe increases. The increase in T corresponds to 
the rise in rigidity of the silo. 

There is one point to be noted in these tables that 
for a =0.5 , p - 0.25 and a =0.5, P = 0.5, OCR and p* is 
more for first case, but e^ is less for second case. 

4.5 MEDira SILO 

4.5.1 Stress Distribution 

Failure in this silo will be same as medium tank, 
that is mode two failtnre. 

The stress distribution and collapse mechanism is 
as shown in Fig. 4,5.1. 

The distribution of nQ is as shown in Fig. 4.5.1(b) . 
There are three zones 

Zone I : Uq = - uSq f or 0 _< e ^ e.| 

Zone II : ng = an^ f or e| < e _< P 

Zone III : Ug = an^ + K;(e-P) n^ for P < e < ^ (4.5.1) 

The boundary and the continuity conditions are 

dm^ j ' 

)^x1 “ = 0 at e = 0 (4.5.2) 
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<^x1 <3.^x2 

m, - = — — ■■■■' — '■ = — ■■ " " e — e 


:1 ~ “x2’ 


de 

dm^2 


de 


1 


dm. 


^x3 


at e = P 


°^2 de de 

Th-e continuity condition at e = 62 depends on vdiether 


(4.5.3) 


(4.5.4) 


If 62 < P 


111^2 “ “®o 


dm 


x2 


— =0 at e = 62 


and 


e? > P 


dm. 


mx3 = - m^ and ' = 0 at e = 62 


(4.5.5) 


= m^ at e = 1 


(4.5.6) 


4,5.2 0 _< e _< (Zone l) 

In this zone, iwhich is from 0 to the 
circumferential stress is negative and reinforcement 
provided is constant and equal to an^. So the equilibrium 
equation (4.3.5) mil be 


d2 

= CO (p*(l-(l+el)"^) + 
de 


First integration gives 


(4.5.7) 


^“x1 


CO (p*(e+T’‘^(l+eT)“*’)+anQe) + A.J (4.5.8) 
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Second integratiocu gjilves 

2 ■ - 2 

es O 

= CG(p* i — + 2-^ln(l+e!P]] --|— ) + L^e+B^ 

( 4 . 5 . 9 ) 

where and atK-e ittonstants of iniesgration. 

From boundary conditions (4.5.2), one gets 

A^ = - CC (4.5.10) 

= 0 (4.5. 11> 

So eq.aations ( 4 . 5 .8 1) and (4.5.9) will be 


de 


= CC (pjjOe+T (l+el)' -S?"" ^)^-anQ^) (4.5.12) 


and 


2 - 2 
e o an^e 

m ^ = GO (p^ ( — + T In (l+&rU- T”"e)+ ----- — ) 


( 4 . 5 . 15 ) 


4.5.3 ®-| 6 ^ P (jtone II) 


In this zone which is from ej| to P, the 
ciroinaferential stress is positive anidfl reinforcement 
provided is const asjjt i„equal to naw e.j to p. So the 

eqnilibri-uia equatiioiai (4.3.5) will be 


2 

dm 


'x2 


= CO (pl'Ctt-d+eO!)"^) -an^jl 


de 


(4.5.14) 


I’irst integration §;|.wes 
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dm^2 

~dS~ "" (pjCe+T" (l+eT)"'^) - ean.^>+-fi 2 <4li,5.15) 

Second integration gives 

e ^ C3tii_ gs ^ 

”>52 = <30 (p« (— + a?-2ln(n.el!):) 

5(44.5.16) 

where and B 2 are constants of int^grraulLon, 

Using the continuity conditions (4-.5.3) , oa.e gete 

Ag = OC (<ie^(n^+H„)-pjT-1j 
OGae^ _ 

®2 = i K ■'■ C)..5.18) 

Substituting these values in equations C: 4 - 5 - 1 : 5 ) aanW ( 4 . 5. 16 ), 
one gets 


de 


OC (pj (e+T-''(l+eT)-''-a:-'')-M^e«e, (ii.V«o» 

( 4..5.19) 


and 


m 


■x2 


SC (p* — + 1 ? ^ln.(l+eT) ’-fl ^ 05 ) 

" ) (M45. 20) 


4.5.4 P < e < 1 (Zone III) 

In this zone, which is from p to i, ties ciocciiiiniferential 
stress IS positive and reinforcement pr-o^i^e.d is nhrearly 
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varying from an^ to from P to 1 . So the equilibrium 
equation ( 4 - 3 . 5 ) vill reduces to 
^2 

= CO (pj(l-(l+eT)“^)-an^-KnQ(e-P)) (4.5.21) 

de'^ 

By first integration 


dm^3 


- OC (p*(e+T“^ (l+eT)"')-anQe - KnQ(^ - Pe)) + 

(4.5.22) 

By second integration 


.-1 


2 3 2 

i»x 3 = OC (p»(^ +0?-‘=ln(1+eT)) - - Eno(h - ^)) 


-2 


+ A,e + 

3 3 


6 2 
(4.5.23) 


where A^ and B^ are constants of integration, 
From the continuity condition ( 4 . 5 , 4 ) , 


A 3 = ac (ae ;(nQ+no) - - P* T"'' ) (4.5.24) 


and 


, KnoP- 
®3 “ ^ "T^ 


ae 


1 


2 


(4.5.25) 


Substitutiiag these values in equations (4.5.22) and 
( 4 . 5 . 23 ), one gets 


'^5 

de 


- CO (p*(e+T""^ ( 1+eT )"" ’-T”^ )-nQae-j+ae^ (nj^+n^) 
■ n„K 


.-1 „,-1 


(e-P)")) 


(4.5.26) 


2 



13 



By using equations (4.5.15), (4.5.20) and (4.5.27), the 
bending moment distribution will be as shown in 
Fig. 4.5, 1(a). 
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4.5.5 Critical Condition 


As tlie value of CO will increase tlie positive 
moment above tbe negative b-inge will increase, and it 
will reach to positive maximum moment niQfo r certain value 
of CC eq[ual to COR at e = e^. This is clear from equation 
( 4 . 5 . 13 ). Ihe value of CCR will also depend on vdiether 
> P . e^ is confined upto 0,1 L of the upper 
portion, so there is no point in considering P < e^. 


So at e = e^. 


dm. 


and 

X o de 


= 0 


( 4 . 5 . 33 ) 


when 


63 < ? 


•^o 


anQ(l+e2l) (t+e-^T) 


and 


2m. 


COR = 


( 4 . 5 . 34 ) 


a3id v^en e ^ > P 


(p*(e^+2I”%n( 1+e^T)-2T””* e^)-anQe|+ (uq+Eq ) ( 2ae^ e2-2i,j ) ) 


n = 


(an^ + n^K/g (e2+e3).-nQKP}(1+e2l)(l+e^l) 


( e ^+ e ^^ + a? 6265) 


( 4 . 5 * 3 ^) 
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and 


6m, 


GGR = 


(3p^(e|+2T"^n(l+e^T)-2T'’^ej) ^5n^ae^ - 

+ aCn^+S^) (6e^e^-3e^)-noK(e3-P)^) (4.5.37) 

4.5.6 Special Case 

let , = £o = = 5o = 1 

Then equation (4.5.28) reduces to 

^ 6/00 f 3a-6a(2ei-6f)tK (1-^) ^ . 

° 3(l+2T"%.n(l4-!r)-2T"'‘) 

When 62 ^ P , equations (4.5.29) and (4.5.30) will be 
ae2 2ae-| . 




(ep + T"”*{l+epT)"’’-l) 


(4.5.39) 


GG = 


(ae|-4ae^e2+2ae^-p*(e2+2T"'%.n(l+&2^)-2T”^e . ) ) 


(4.5.40) 

when 62 2 P » equations (4.5.31) and (4.5.32) will be 


P* = 


2ae2 + K(e2-P ) ^40^6^ 


.2(65+1“'’ ( 1+621)" ‘ - T“‘ ) 


-1 


(4.5.41) 


and 


GO = 


(3ae^-f 6aef +5(e2-P)^-1 2ae^ 62 -p*( 3e|+6T“^(e2T+1 ) 

--6T”‘’ep)) (4.5.42) 
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vdien _< p , equations (4.5.34) and (4.5.55) will be 

= a(l+e2^) U+e^'^) / (4.5.43) 


and 


OCR = 


p*(e^+ 2 T”^ln(l+e 52 !)- 2 I"''e^)-anQe|+ 2 ( 2 ae^e 2 -ae^) 

(4.5.44) 

when > P , equations (4.5.36) and (4,5.37) will be 


PS = 


((a+E:/2(e2+e^)-E:p)(l+e2g)(l+e3!r)) ^ 

2! (e2 6^ ®2^3^ 


and 


OCR = 


P*(3e|+62!”%n(l+e5!r)-6T"'‘e^.3a€>| 

-K(e5-P)^+6a(3e^e3-e^)) (4.5.46) 


when a = p =1 


p* = (6/GG+3-6e^(2-e^))/ 3(l+2T“^ln(l+T)-2T”'’) 

(4.5.47) 

Pq = (e2-2e-,)/(e2+^C“l(l+e22l)“*‘-l)^^^^^^^^^^^^^^^^^^^^^ (4.5,48) 

CG . = 2/(e|-2(2e^e2-e^>^)^(e2+2T’^(l+e2^)-2T"‘^e2)) 


(5.4.49) 



77 


Pj = (l+e2T)(l+e^T)/(e2T+e^T+e2e^T^) (4.5.50) 

OOR = 2/(p*(e|+2T’'^ln(l+e3T)-2T“''e^)) (4.5.51) 

4.5.7 Critical Yalues 

Tlie distribution of ng is dependent upon the values 
of a and P. For a given set of a and P values at certain 1 
the critical shell parameter value OCR, location of hinges 
e 2 and e^, the range of hoop compression from 0 to e^ , and 
collapse pressure p* are shown in Table 4.4 for m^ = 51^= n^ 
= Uq = 1 and T=1. 

TABLE 4.4 : COR E OR SECOUL MOLE OF FAILURE OR SILOS AT 



T = 

1 ANL m^ 

II 

II 

“o = 

> = 1 


a 

P 

. OCR 

®i 

®2 

®3 

PS 

0.25 

0.25 


■ 




0.50 

0.25 

- 

- 

- 

- 

- 

0.50 

0,50 

- 

- 

- 

- 

- 

0.75 

0.50 

- 

- 

- 

- 

- 

1.00 

1.00 

197.2 

0.061 

0*64 

0.214 

2*006 
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From Table 4.4 it is clear that the OCR for second 
mode, for a an.d p less than, one ,is more than 400 . So the 
range of second mode is large and as the value of T 
increases this range will also increase. 

4.6 DEEP SILO 

4.6,1 Stress Distribution 


Failure in this mode is known as mode three failure 
and also known as partial failiire or partial collapse because 
in this failure upper portion remains intant as in the third 
mode of failure in tank. The stress distribution and 
collapse mechanism are as shown in Fig. 4.6.1. The 
distribution of ng is as shown in Fig, 4.6.1(b), This is 
given by the equations. 


Zone I : ^9 = ^ ® ^ P 

Zone II : ng = anQ+K(e-P)nQ for P ^ e ^ 1 (4.6.1) 

The boundary and the continuity conditions are 

°^3t1 = e = e4.for < P 

%2 ~ “0 » =0 at e= f or e^ > p ( 4 . 6 . 2 ) 



= ®x2# 


^1 

de 


at e = P for e^X p ( 4 . 6 . 3 ) 


^%2 

de 
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The continuity con'^itions at e = depends on blether 
®5 ^ P 


5 'or ©5 ^ P 

dm 

- “®o » ° at e = 63 

and for 2 P 

_ ^x 2 

“x2 = - “0 - = 0 at e = 65 

®^d m ^2 = at e = 1 


(4.5.4) 

(4.6.5) 


4.6.2 < e < p (Zone l) 

In this zone, circumferential stress is positive and 
reinforcement is constant throughout, from e^ to p, is an^. 
So equation ( 4 . 3 . 5 ) reduces to 

d^m^-i 

= GO (p*(l-(l+eT) ”*) - an^) (4.6.6) 

First integration gives 

dm^l 

= CC (p* (e+T*“''(l+eT)“'')-anQe) + (4.6.7) 

Second integration gives 

2 2 
nixl ~ (Pg(~ + T~^(l-heT))- - )+ 

(4.6.8) 
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where and are constants of integration. 


Prom boimdany conditions (4.6.2) 

= - GC (pj(e^+T”''(l+e4T)’''')-anQe^) 


B, = 00 (p* (~^ -T ln(l+e.T)+[I!" e.d+e/)" )- — ) 


'4 


o ^ 2 


+ 21 ^ (4.6.9) 


Substituting these values in equations (4.6,7) and ( 4 . 6 . 8 ), 
one gets 


dm. 


■x1 


= 0G(p*(e+T“^ (l+e!]?)”'-e^-T"' (l+e^T)"''-anQe+anQe^) 

( 4 . 6 . 10 ) 


.-1 ^-1 


.-1 


de 


and 


m. 


= CC (pg((e-e.^)V2 +T'^(ln(1+eT)-ln(1+e^T)). 




'^/ / ‘ 


4.6.3 P < e < 1 (Zone II) 


In this zone the circumferential stress is positive 
from p to 1 , and the reinforcement provided in this zone is 
linearly varying from an^ to n^. So the equation (4.3.5) 
will be 


de^ 


00 (p*(l-(l+el)‘"^)-anQ-Enp(e-p ) ) 


( 4 . 6 . 12 ) 
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Integrating once 

diii„ p 11 ? V 

-"^ = 00 (p*(e+l" (1+eT)” )-anQe-ln^(eV2-^e))+A5 


de 


( 4 . 6 . 15 ) 


Integrating once again 


™x2 = (pJ(eV2+T"^ln(l+eai))- an^eVs -Kn^CeVS-PeVs) ) 


+ A^e + 


(4.6.14) 


where A^ and are constants of integration. 


Brom the continuity conditions (4.6.3) 

A5 = CO (anQe^-p*(e^+T-^(H-e4^r)”'‘)-KnQpV2) (4.6.15) 

and 

B 5 = CC(p*(e|/2-5'*“^ln(l+e^T)+e^T“'’(l+e^T)“^) 

-anQe|/2 + KnQpV6) + (4. 6 . 16) 

Substituting these values in equation (4.6. 1 3) and (4.6. I4) , 
one gets 


d^ 

de 


CG(p*(e+l"‘''C(1+eI r2 (1+e^I )“Ve^) 
-anQ(e-e^)- 111^/2 (e-p )^) 


and 


(4.6.17) 


111^2= GG(p*((e-e^)^/2+ai”^(ln(.1+eT)-ln(l+e^T))- 

-.T”^(l+e^T) ^ (e-e^))-JDcnQ/2(e-e^)^-Enjj/6(e-P )^)+mQ 

( 4 . 6 . 18 ) 
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using the boundary condition ( 4 . 6 . 5 ) in equation ( 4 . 6 . 18) 
3ano(l-e.)^+Kno(l-p)^ 

° 5 ( ( 1 -e^)^+2T“^ (ln( 1+!E )-ln( l+e^^T ) )-2T"‘‘ ( 1+e^T )’"^ ( 1 -e^) ) 

(4.6.19) 

4.6.4 < P (Case l) 


For ^ P » using the continuity conditions (4,6,4) 
in equations (4.6,11) 


p* = 


ano(e^-e^) 


(4.6.20) 


"" (e^-e^) + T”'‘((l+e^T)"''-(l+e,a})"'') 


and 


GO = 


2(mQ-t^o) 


an^ (e^-e^) ^-p* ( (e5-e^)^+2T*"^ (ind+e^T )-ln( 1+e^T ) ) 


— 2T (l+e^T) (e^— e^)) (4.6.21) 

For e^2 P , using the continuity conditions (4.6.4) in 
equations (4.6, t7) and (4. 6. 18), one gets 




and 


GC = 


- ^ (^.6.22) 
2((e5-e^) + !T”'‘((l+e5T)"'‘-(l+e^T)-'') ) 

6(mQ+5o) , 

3ano(e5-e^)^+KnQ(e5-P)^^p*(3(e5-e,)^+6!E’'^(ln(l+e5'3}) 


-ln(l+e^T))-T“'‘(l+e^a))^'‘(e5-e^)) (4.6.23) 
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4.6.5 P < (Case II) 

lor this case equilibrium equation (4.3.5) is 
a m^-r 2 

— = 00 (p*(l-(H-eT) -an -Kn^(e-P)) ( 4 . 6 . 24 ) 

deZ ° 

Integrating once 

— = 00 (p*(e+T'’^(l+e'E)“^)-an-e-^hio(eV2-Pe)) + Ag 

( 4 . 6 . 25 ) 

Integrating once again 

2 

mx 3 = 00 (p*(eV2+T’“'^n(l+el) )- - 5io (eV6-PeV2) ) 

+Age. 4- Bg (4.6.26) 

Where Ag and Bg are constants of integration. 

Using continuity conditions ( 4 . 6 . 2 ) , one gets 

■^6 ~ “0C(p*(e^+T (l+e^T) — anQe^-KnQ(e^/2-Pe^)) 

( 4 . 6 . 27 ) 

and 

Bg = m^ - CC(p*o(T“^ln(l+e^T)-e|/2-I''''e^(l+e^T)”‘’ ) 

2 

+ ~ 2 §e^/2) (4.6.28) 

Substituting these values in equation ( 4 . 6 , 25 ) and (4.6.26), 
one gets 
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am ^3 

Te 


and 


= OC(p*(e+a?“'' (l+e!E)"''-e^-T"'^ (l+e^T)”’’)-anQ(e-e^) 


-En^CeVs - Pe - e|/2 + pe^^)) 


o- 

(4.6.29) 


m 


(e-e r 

x3 = CG(p*( ■■ ■ — + T"^(ln(l+eT)-ln(l+e4T)) 


(4.6.30) 


-T"' (l+e^T)" (e-e^)) -anQ/2(e-e^) -En^/e 

(e^-3Pe^+3e e^-6pee^+2e^-3Pe^)) + mo 

From boundary condition (4.6.5) 

3 ano(l-e^)^+Eno(l- 3 P+ 6 pe^- 3 e^- 3 Pe^ + 2e^ ) 

3 ( (1 -e^) ^+22“^ (ln( 1+T )-ln( ) )-2T"'‘ ( 1+e^F)"'’ ( 1 ..e^) ) 

(4.6.31) 

From the cantinnity conditions (4.6.4) 

_ 2gno(eg-e^)+Enp(e|-2e^P-e4+2peP ,(4..6.32) 

° ” (e,+T-''((1+e.T)-''-(l+e,Tr'') - e. ) 


P# 

■^o 


GO = 


3^-^ y l-rc^J 

6(mo+mo) 


3ano(e5-e4)^+Eno(e|-3Pe|+3e5e|-6pe^e^+2e|-3Pe^ ) 

- 3 p*( (e5-e^)^+2T“^(ln(l+e50: )_ln(l+e4T ))-'!“'' (l+e^F)“^ (e^-e^)) 


(4.6.35) 


4.6.6 Special Case 

let mo = 5o = no = So = 1 
For e^ X P, equation (4.6.19) will be 
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)* = 


3 a(l-e^)^ + £( 1 -P)^ 


° 3((l-e^)^+2T"^(ln(l+T)-ln(l+e4^T))-T"'’ (1+64^3?)”^ (1-e^)) 

( 4 . 6 . 34 -) 

when ^ p , equations (4. 6. 20) and (4.6,21) will be 

p» = (+.6.55) 


° (e 5 -e 4 )+!D"''((l+e 5 T) '^-(l+e^T)"'* ) 


and 


CO = 


a 


(e5-e4)^-p*((e5-e4)^+20:”^(ln(l+e5T)-ln(l+e4T)) 


- 2 T“''(l+e 4 l)"''(e 5 -e 4 )) 


(4-6.36) 


when ©5 2 P » equations (4.6.22) and (4.6.23) will be 

2 a(e 5 -e 4 ) + K(e^-P)^ 

'b ^ ^ 


Pf = 


+ 2 T”''((l+ep;T)'''^-(l+e,T)“'' ) 


(4.6.37) 


and 


12 


GO = 


( 3a ( 65-64) ^+K ( e^-p ) ( 3 ( e^-e^) ^+ 61 ”^ (in {A+e^T ) 


5 -= 4 > 


-ln(l 4 -e 4 l))-T""^ (1+64!) ^ (6^-64) )) 


( 4 . 6 . 38 ) 


vdien 64 2 P t the equations (4.6.31), ( 4 . 6 ^ 32 ) and 
(4.6.33), mil be 


P* = 


^(1-64)^ + K(l- 3 P+ 6 pe 4 - 3 e|- 3 pe|+ 2e|) : 

° 3 ((l-e 4 )^+ 2 T“^(ln(l+T)-ln(l+e 4 l) )-2f”‘* (1+641)“'* (1-64)) 

( 4 . 6 . 39 ) 
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Three graphs have been drawn between GO and p* for 
different values of T, a and P, The S'ig* 4*7.1 has been 
drawn for a = p = 1 and for different T values. It is 
clear from this graph that as the value of T increases the 
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C =4^/tR 
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pressure' ■ d i s-t r' j'b u t i ■■ 
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range of CG for first mode increases. When the critical 
points of all the 1 are joined for a = p = 1 , one gets a 
curve which separates the mode I from mode II is 

158.2258 2 

p* = - 6.2098 + 0.1237 COR + - 0.000? OCR 

° CCR 

(4.7.1) 

In Rig, 4,7.2, curve a = p = T = 1 intersects curve 
oc = 0,75, P = 0,5 and 1 = 1 , at GO = 208. This is because 
the modes of failure for these two cuiives are different. 

The former is in the third mode whether later is in the 
second mode. As e 2 is always less then p, the curves are 
smooth and changes in slope are only at the critical point 
e^ is always confined to 0,075 of the upper portion for 
GG < 400 and T = 1. 

In Rig, 4. 7. 5 for a =*■ p = 0.25 and T = 2,5 there 
is only one mode of failure for CG < 4 OO. Rrom Rig. 4.7.1 
it is clear that as the value of T increases the value of 


e 2 decreases 



5. ANALYSIS OP REIKPORGED GCSfCRETE SJCLOS BASED ON 
THE JANSSEN'S PRESSURE DISTRIBUTION THEORY 

5.1 PRESSURE DISTRIBUTION 

TNe material filled in cylindrical silos, in vdiich 

emptying Hole situated at tlie centre of the bottom, is 

held in equilibrium by the reaction from the sloping bottom 

and the frictional force on the vertical wall, as shown 
( 1 ■5) 

in Pig. 5.1.1. . Equations are given below. 

The vertical static unit pressure at depth X below 
the surface of the stored materia;V is 

<1 = ,-r-^ ( 1 - exp ( -p-iE 2/an ) ( 5 . 1 . 1 ) 

M- k * 

The lateral static unit pressure at depth X is 
Y R 

P T" (1 - e3:p(-4i’kX/R^)) (5.1.2) 

Yertical frictional force per unit width of wall perimeter 
on wall above depth X 

Y = (Y X - 0.8 q ) R.J (5.1.5) 

where 

p' = tan coefficient of friction between stored 
material and silo walls ,and 
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1 - sin 0 

]c = — , 

1 sin / * 

vfcere 

0-* = angle of friction between material stored and 

silo wall 

0 = angle of internal friction of stored material 

= hydraulic radius 
Y = unit weight of stored material 

Static unit pressure normal to a surface inclined 

at an angle 6 to the horizontal at depth X below surface 

( 12 ) 

of stored material is^ 

2 2 

= P sin 9 + q. cos 6 (5.1.4) 

5.2 EQUILimrJl! EQUATIONS 


The equilibrium equations for a silo in non- 
dimensional form are given by the equations (2,5.5) and 
(2, 3. 6)* Unlike water tanks there is a vertical pressure 
along the walls besides the normal horizontal. The 
equilibrium equation for the normal pressure is given by 
tR 

+ ng = p* (5.2.1) 


41.' 


^ % 
de^ 


where 



2 YRi , ^ • 

= - — 1™. (l-e 2 pC-kLix»X/(R.l))) 

p’oi,t 

= (l-expC-rTe)) 


( 5 . 2 . 2 ) 
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A shallow silo is also known as short silo. The 
failure of the silo at collapse is known as more one failnre 
as in tank. The distribution of the stress resultants and 
the collapse me chani an are shown in Pig, 5.5.1. The 
distribution of Ug is as shown in Pig. 5.3.1(c), This is 
given by the equations 


Zone I : ii0 = “Hq for 0 _< e _< p 

Zone II : jiq = an^ + Kn^Ce-P) for P < e X 1 (5.3>1) 

vUiexe 

1 -« 

^ “ T^r 

The boundary and the continuity conditions are 

d-TTLy. 

° ^ ^ ® ° (5.5.2) 
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= m. 


and 


"x1 


“x1 


‘x2 ’ 


de 


dm 


— — at e = p 

de 


= m at e = 1 


5.3.2 0 < e < p (Zone l) 

In this zone, which is from 0 to P, 
reinforcement is constant and equal to an^. 
equation (5.2.5) reduces to 

dm_., 

— ^ = CC (p*(l-exp(-Te)) 
de^ 

Integrating once 

dTTl j^ ^ ^ 

= CC (p*(e+T~‘exp(-Te) ) -an^e) + i 

de ^ 


Integrating once again 

2 m — 2. 


an^e 


114.^ = CC(p*(eV2-T"'^exp(-Te)) - 


where and are constants of integration, 

Using boundary condition (5'. 3.2) , one gets 
A,, = - CC p* 1“'' 

and B,j = CC P^”^ 

So equations (5.3.6) and (5.3.7) will be 


(5.3.3) 

(5.3.4) 

the circumferential 
So the equilibrium 

(5.3.5) 

i.-| (5.3.6) 

)+A^e+3',, 

(5.3.7) 

(5.3.8) 

(5.3.9) 
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— ■ = CC (p*(e+T”"*exp(-Te) -5}“"* )-arL^e) 


(5.3.10) 


and 


m 


xi = OG (p*(eV2-T~^exp(-!re)-!P"‘’e+T’^)-anQeV2) 

(5.3.11) 


5.5.3 P < e ^ 1 (Zone II) 


In this '7;one, vhich is from p to 1 , the circumf erential 
reinforcement is varying linearly from an^ to n^. So the 
equilibriiim eq.uation (5.2.5) reduces to 


d^m^2 

^^7“ =00 ( ( 1 -exp ( -T e ) ) ( e -P ) ) 

Pirst integration gives 


(5.3.12) 


dm^2 -t 

= CC(p*(e+l“ exp(-Ie))-anQe-KnQ(eV2-Pe)) + A 2 

(5.3.13) 

Second integration gives 

m^2 = (p*(e^/2-T“^exp(-Te))-anoeV2-^o(eV6-Pe^/2)) 

+ A2e + B2 (5.3.14) 

where A 2 and B 2 are constants of integration. 

The continuity conditions (5.3.5) gives 


A 2 = - CC (T"''p* + EnQp^/2) 
and B 2 = CC (T’^p* + EiioP ^/6 ) 


( 5 . 5 . 15 ) 

(5.3.16) 
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So e(iuations (5.3.13) stad (5.3.14) will be 


diiL_p -j -I p 

= 00(p*(e+T~ exp(-!re)-T'" )-anj^e-Kn^/2(e-P ) ) 


de 


aind 


m, 


(5.3.17) 

:^2 = CC (p*(eV2-T~^exp(-Te)-T“”'e+T''^)-anQeV2-KnQ(e-§)V6) 

(5.3.18) 


Hence from boundary condition (5.3.4) 
em^/GC + 3anQ+Kn (l-p)^ 

° 3(l+2!P"^(l-exp(-T))-2T~-’) 


(5.3.19) 


By using tHe equations (5.3.11) and (5.3.18), the bending 
moment distribution will be as shown in ?ig. 5.3.1(b), 

5.3.4 Critical Condition 


As the L/^l ratio of the shell increases, the value 
of GG as well as 1 increases. If l/R ratio renains same, 
but L and R increase , the value of CG increases but I 
remains constant. As 00 increases, p^ decreases, so m^-| 
and m ^2 will also decrease and it will reaeh the maximimi 
negative value for certain 00 = OCR and certain e = ©q* 

So for e = e^ 

“x ~ ” ®o dm^-Zde = 0 (5.3.20) 

when e^ X p 
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(5.3.21) 


an^e J-p* ( 63+ 2T"^ ( 1 -exp ( -T e^ ) ) -2T" ' e^ ) 


(5.3.22) 


when P ^ ©o ^ ^ 


Vt = 


2an^e^+-Sji^ieo-^y 


° 2(e„+!P"^exp(-Te„)-T“'‘) 


(5.3.23) 


OCR = 


3anQe^+Kno ( e^-P ) ^-p* ( 3eo+6I"^ ( 1 -exp ( -Te 0 ) ) -60}"'’ ) 


(5.3.24) 


5.3.5 Special Case 


let = mo = = 1 


Olhen 


6/CC+3a+K(l-P )^ 


P* Z= 

° 3(l+2T"^(l-exp(-a})-2T”'') 


when Bq ^ P 


= aeo/(eo+T” exp(-Ten)~T” ) 


COR = 


(5.3.25) 


(5.3.26) 

(5.3.27) 


aeo-p*(eo+20}“^(l-exp(-0}eo))-2a!”''eo ) 


when P ^ Sq ^ 1 

2 aeo + E (e„-P )2 

^TT — ■. Ill I ■■■■«■■■ I I nail 

° 2(e^+a}"^exp(-Teo)-0!”^ ) 


(5.3.28) 
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and 


OCR = 


3aeQ+K(eQ-P )^-p* C 3 e^+ 6 T “^ ( 1 -exp (-TeQ ) )- 6 T'"'’ ) 

(5.3.29) 

when a = p = 1 

1 + 2/CC 


PS 


(1 + 2T"^(l-exp(-2))-2!I!“‘) 


-1 


(5.3.30) 


PS = 


(Oq + exp(-a?eQ)-T“^ ) 


(5.3.31) 


and 


COR = 


( e^-p* ( ( 1 -exp (-Te^ ) )- 2 T“ ' e^ ) 


-1 


(5.3.32) 


5.3.6 Critical Values 

The distribution of n 0 is dependent upon the ■values 
of a and P. Ror a given set of a and P values and for pre- 
cribed T, the critical shell parameter value OCR, the 
location of hinge e^^ and collapse pressure p* are sho'wn 
in Tables 5.1, 5.2 and 5.3 for = nQ=nQ =1 and for 

T eq,ual to 1, 2.5 and 4 respectively. 
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TASLE 5.1 

: OCR POR FIRST MODS OF 

AT T = 1 AND 111^ = 5^: 

FAILURE 

= no = 1 

OF SILOS 

a 

P 

GC5R 

®o 

PS 

0.25 

0.25 

169.00 

0.4565 

1.5252 

0:50 

0.25 

49.20 

0.5497 

2.4014 

0,50 

0.50 

45.80 

0.5358 

2.2178 

0,15 

0.50 

27.75 

0.5628 

3.1965 

1.00 

1.00 

19.90 

0.5719 

4.1628 


TABLE 5.2 

: OCR FOR FIRST MODS OF FAILURB 

T = 2.5 AKD fflo = Sq = n^ = 1 

OF SILOS AT 

a 


GGR 


Po 

0.25 

0.25 

266.0 

0.335 

0.8064 

0.50 

0. 25 

67.5 

0.498 

1.2628 

0.50 

0. 50 

59.2 

0. 499 

1.1653 

1.00 

1 . 00 

24.95 

0.5479 

2.1874 
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TABLE 5.3 : OCR EOR FIRST MODE OF FAILURE OF SILOS AT 



T = 

4 and m^ 

11 

o 

II 

11 


a 

P 

GGR 

®0 


.0.25 

0.25 

382.0 

0.276 

0.6357 

0.50 

0.25 

90.2 

0.447 

0.9887 

0.50 

0.50 

75.0 

0. 468 

0.9125 

0.75 

0.50 

43.5 

0.507 

1.3125 

1.00 

1.00 

30.5 

0.5259 

1.7129 


From these tables it is clear that as the value of T 
increases the range of GGR for first mode increases. T 
increases means L/k increases so the rigidity increases. 

5.4 MEDIUM SILO 

5.4.1 Stress Distribution 

Failure in this silo will be of mode two type as 
in medium tank. The stress distribution and collapse 
mechanism are as shown in Fig. 5.4.1. The distribution 
of n^ is as shown in Fig. 5.4.1(b). There are three 
zones: 

Zone I : n^ = - an^ for 0 < e ^ e^ 

Zoiae II : n.Q = an^ for e-j ^ e ^ p 

Zone III : ■ n0 = an^+KCe-p) for p _< e X 1 


(5.4.1) 
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The boundary and the continuity conditions ane 


‘x1 

dnu.^ 

= 0 , — = 
de 

0 at e = 0 

( 5 . 4 . 2 ) 



- 1^+ 0 - 0 

( 5 . 4 . 3 ) 

‘x1 

^2 ’ (3^0 

de ^ 

x2 

dm ^2 

”^3 ’ de 

dm„5; 

_ — ^ at e - P 
de 

( 5 . 4 . 4 ) 


The continuity condition at e = 02 depends on vdiether 



If 62 < P 

diiijj2 

m_p = hE^ and 0 at e = e^ 

AC u de 

and for 62 > P 

in ^3 = -mQ and — = 0 at e = e 2 (5.4.5) 

and 

m ^3 = at e = 1 (5.4.6) 

5 . 4.2 0 ^ e ^ e^ (Zone I ) 


In this zone, which is from 0 to e^ , the circumf erential 
stress is negative and reinforcement provided is constant and 
equal to ^5^. So the equilibrium equation (5.2,5) will be 


de2 


GO (p*(l-exp(-Te))+ an^') 


(5.4.7) 


Integrating once 
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dm^l 


= GO (p*(e+T‘"”*exp(-Te) ) + ocnge) + (5.4.8) 


Second integration gives 

m ^1 = GO exp(-Te)) + anQe^/2) +A^e+p^ 


where A^ and are constants of integration, 
Brom the boundary conditions (5.4.3), one gets 
A^ = - CG p* 

= GC p^"^ 

So eq,uation (5.4.8) and (5.4.9) reduce to 


(5.4.9) 


(5.4.10) 

(5.4.11) 


1 TTl .j A A 

—22- = GG(p*(e+T” exp(-Te) -T"') + aS^e) 
de ° 


(5.4.12) 


= GC(p*(eV2+T“^(l-exp(-Te))-T”‘’e)+ anoeV2) 


(5,4.13) 


5.4.3 i e < P (Zone II) 


In this zone, which is from e^ to p, the circumferential 
stress is positive and reinforcement provided is constant 
equal to an^ from e^ to p. So the equilibrium equation 
(5. 2.5) mil be 


= GG (p*(l-exp(-Te)) - an^) 


(5.4.14) 



107 


I'irst integration gives 


dm, 


de 


■ - GG(p*(*e+!D“^exp(-Te))-an^e) + A 2 


(5.4.15) 


Integrating once again 


m^2 “ GG(p* (eV2-^”^®^P(-^e))- an^e^/^) + k2^+'B2 

(5.4.16) 

where A 2 and B 2 are constants of integration. 

Using the continuity conditions (5.4.3), one gets 

A 2 = CG (-p^"’Vae^(n^+no)) (5.4.17) 

and 

B 2 = GG(pgT"^-aef/2 (n^+Ho)) (5*4.18) 

Substituting these values in equations (5.4.15) and (5.4.16) 


dm^2 1 1 ' 

— = GO (p*(e+T’’ exp(-Te) -T“ )-.anQe+ae-|(n^+nQ)) 

( 5 . 4 , 19 ) 

and 

®x2 = 0G(p*(e^/2-.l“^exp(-Te) -T"''e + T’^) 

-anQeV2+(aee-,-ae^/2) (Uq+Uq)) (5.4.20) 

5 . 4.4 P < e < 1 (Zone III) 

In this zone, which is from P to 1, the 
circumferential stress is positive and reinforcement provided 
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is linearly varying from an^ "to P to 1, So the 

equilibiiiaii equation (5.3.5) will be 


de2 


GC(p*(l-exp(-le)) - anQ-Kn^jCe-P )) (5.4.21) 


First integration gives 


de 


- CO (p*(e+T“^Sxp(-Te) ) -anQe-En(^(e^/2-Pe) )+A^ 

(5,4.22) 

Second integration gives 


“x3 ~ exp(-Te))-anQeV2-KnQ(eV6-Pe^/2) ) 

+ A^e + (5.4.23) 

where A^ and are constants of integration. 

Using the continuity conditions (5.4.4). 

Aj = - 00(1“'' p* + KjioPV2«e^(n^+5„)) (5.4.24) 

and 

B 3 = GO (p^‘^+ Kn^p V6-aef/2(nQ+nQ)) (5.4.25) 

Substituting these values in equations (5.4.22) and 

(5.4.23). 

1 ^ o 

— = CG(p* (e+T" exp(-a}e)-l"'^)=anQe-KnQ/2(e-P) 

^ 0 

+ae^(n^+E^)) 


(5.4.26) 
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= CC(p*(e^/2-T”^exp(-Te)-T"e+T''^)-anQeV2 

KnQ/6(e-P)^+(nQ+nQ)a(e-je-e|/2)) (5.4.2?) 

Using the boundary condition (5.4.6) in the equation 
(5.4.2?) 

6mp/GC+5ano+£np(l-P)^-5a(nQ+no)(2e.^-ef ) ^^ ^ 

° 3(l-2T“^exp (-T) -2T“V2T''^ ) 

From the first condition of (5.4.5) , one gets 

Q -^ o®2 ~ *^^1 ^^ o "^ ^ o ^ ^ ( 5 . 4 . 29 ) 


P * = 


(e^ + T”"* exp (-Tep) - T”' ) 


-1 


and 


25 , 


GC = 


ry 

(an^el" (n^+n^ )a ( 2e^ e2-ep -p* ( e2-2T”'^exp ( -162 ) 


-21 "* 62 + 21 ”'^) 


- 2 ' 


( 5 . 4 . 30 ) 


From second condition of (5.4.5) 

2 — 

2an e 2 +Eno(e 2 -P) -ae^(nQ+nQ) 

■ n * — ,,,, — i , i „ .,, 


and 


GG = 


2 (e 2 + 2 !”' exp (- 162 ) - ) 

6 mo 


-P*(3e2-6T"^exp(rTe^T“”'e2+6T"^)) 


( 5 . 4 . 31 ) 


( 3anQe2+£nj^ ( e 2 “P ) - 3 a ) ( 2e^ ® 2"'®1 ^ 


( 5 . 4 . 32 ) 


By using equations (5.4.13), (5.4.20) and (5.4.23), the 
bending moment distribution will be as shown in Fig.5.4.1 (a) 
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5.4.5 Critical Condition 


As the value of CC increases the positive moment 
above the negative hinge will increase and it will reach 
the positive maximum moment m^ for certain value of 
CO equal to OCR at e = e^^. The value of CCR will also 
depend on vdiether > P. e^ is confined to if^thin 0, 1 L 

of the upper portion so there is no point in considering 

p < e-i . 


So at e = e- 


“x = “o ‘“'5 37 


(5.4.33) 


when 


< P 


P* = 


anQ(e2-e,) 


( 62 - 6 ^) + T””’ (exp(_Te 2 ) -expC-Te^)) 


(5.4.54) 


CCR ~ 2 P 1 2 

(p*(e^+2T"^(l-exp(-Te5))-2T" e^) -an^e^ 

+a(n^+no)(2e.,e^-ef)) (5.4.55) 

when P ^ ©3 ^ ©2 

p* (2anQ+Kno ) (e2+e3)--2£noP ) ( 62 - 53 ) (5.4.36) 

° 2((e2-e^)+!E*^exp(-!Pe2)-exp(-Te3)) 

and OCR = 6mo/(P^(5e5+6T‘’^(t-exp(-Te5))-6!D‘’^e5)-5anQe| 

-Kno(e3-p)^+5a(no+EJ,)(2e^e^-e^)) (5.4.57) 
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5.4*6 Special Case 


let fflo = hIq = = Hq = 1 

Hen equation (5.4.28) will be 

6 /CC + 3a + K(l-p)^*' 6 ae^ (2 -e-j ) 

° ( 3 + 6 T"^(l-exp(-l)) - 6 T“^ 

v4ien 62 ^ P , equations (5.4.29) and (5.4.30) will be 


(5.4.38) 


p* = 


a(e2-2e^ ) 


(e2 + (exp(-!Ee2)- I)) 


( 5 . 4 . 39 ) 


and 


OG = 


(ae|-2a(2e^e2-ef)- p*(e|-2!I?”^exp(-Te2) ;) 

( 5 . 4 . 40 ) 

when 62 4 I P , equations (5.4.31) and (5.4.32) will be 


p* = 

■*^0 


and 


2 ae2 + K(e2-P)^ -4ae.^ 

2(e2 + 3)“"* (exp(-Te2) -I)) 

6 


(5.4.41) 


CG = 


(3ae|+K(e2-P)^-6a(2e^e2 -ef)-p*(3e|+6T"^(l-exp(Te2)) 
- 6T"^e2)) ( 5 . 4 . 42 ) 

when _< p , equations (5.4.54) and ( 5 . 4 . 35 ) reduce to 


2 -2 


n = 


a (e2”"®^) 


( 5 . 4 . 45 ) 


i 

((62-65) + 1“ (exp(-Te2) -^expC-Te^)) 
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and 


GCR = 


(p*(e^+2!E"^(l-exp(>Te^‘)-. Te^) )- ae^+2a(2e^e^-e:^) ) 

(5.4.44) 

vdien P _< eq^uations (5.4.36) and (5.4.37) will be 
( 2 a + K ( e 5 + e ,)- 2 Kp )( e ^- e „) 


p* = 

■^o 


-2 ^ 3 ' 


■?,: 


2 ((e 2 --e^) + T”^ (exp(-Te 2 ) -exp(--le^)) 


(5.4.45) 


and 


OCR = - 


( pS ( 3 e 5 + 6 T "^( l - exp (- Te ^) T ' e ^)) 

+ 6 oc ( 2 e ^ ® 3'"®1 ^ ) ( 5 . 4 * 46 ) 

when a = p = 1 

2/GG + 1-2e^ (2-6^ 


P ^- = 


(l+2T“^(l-exp(-T))-2!E''' ) 
( e 2 - 2 e ^) 


-1 


•^o 


GO = 


(e 2 + T“^exp (-162) ) 

2 


,-1 


(5.4.47) 


(5.4.48) 


(e|-2(2e^e2-e^)-p$(e^+2a}"^(l-exp(Te2))-2T""’e2 ) 


p* 


( e ^- e 3 ) 


A 

((e 2 -e^) + 1 (exp(-Te 2 ) -exp(-Te^ )) 


(5.4.49) 

( 5 . 4 . 50 ) 


CGR = 


(p* (e|+2T"2 ( 1 -exp (-Te 3 ) ) - 23 } •'*e^ ) -63+2 ( 2e^ e^-ef ) ) 

( 5 . 4 . 51 ) 



113 


5.4.7 Oritical Values 

The distribution of Uq is dependent upon the values 
of a and p. For a given set of a and P values at certain T 
the critical shell parameter value COR, location of hinges 
62 and e^, the range of hoop compression from 0 to e^ and 
the collapse pressure pg are shown in Tables 5.4 and 5.5 for 
= Sq = n^ = n^ = 1 and T equal to 1 and 2.5 respectively. 


TABLE 

5.4 : OCR FOR 

SECOND MODE OF 

FAILURE 

OF SILOS 

AT 


T = 1 AND m^ 

II 

II 


1 


a 

P 

OCR 

®1 

e2 

®5 

pS 

0.25 

0.25 

— 

— 

— 

* 

-- 

0.50 

0 . 25 

- 

- 

- 

- 

- 

0.50 

0.50 

- 

- 

- 

- 

- 

0.75 

0.5 

299.0 

0.065 

0.595 

0.212 

2.26 

1.00 

1.0 

119.5 

0.075 

0,750 

0.315 

2.70 


TABLE 5.5 : CCR FOR SECORD MODE OP FAILURE CF SILOS AT 



T=2. 

5 A3JD 


If 

II 

So = 1 


a 

P 

OCR 

«1 

®2 

63 

PS 

0,25 

0.25 

— 

y 1 ■ 1 1 r 1 . 


— 

— 

0,50 

0.25 

- 

- 

- 

- 

rr 

0^25 

0,50 

• - 

- 

- 

- 

- 

0.75 

0,50 

92 

0.041 

0.55 

0.202 

1,36 

1.00 

1.00 

275 

0.060 

0.525 

0.288 

1,62 


As T increases, tlie zone for second mode for 
a * 0.75 and P = 0.5 decreases. 

5.5 LISP SILO 

5.5.1 Stress Distribirtion 

Failnre in this mode is known as mode three failtire 
and also known as ^partial failure’ or ’partial collapse*, 
same as in the case of tank. 

The stress distribution and collapse mechaniam 
are as shown in Fig. 5,5.1. 

The Idi^ribution of n0 is as shown in Fig, 5. 5.1(b) 
This is given by the equafriohs 
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Zone I : Hg = an^ for < ® < P 

Zone II : ng = an^ + K (e-p) nQ for P ^ e ^ 1 (5.5.1) 



In th.is zone» circumferential stress is positive 
and reinforcement in constant throu^out, from e^ to 
is an^ so equation (5-2.5) reduces to 
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a\^/cle^ = CC(p*(l-exp(-Te))-ano) (5.5,6) 

Integrating once 

dm^-j/de = GC(p*(e+T‘"^exp(~Te))-anQe) + A-j (5.5.7) 

Integrating once again 

= CC (p* (6^2-1"^ exp(-Ie))-anQeV2) + A^e+B^ 

(5,5.8) 

Tf^ere A-j and are constants of integration. 

Prom the boundary conditions (5.5.2) 

A-i = - 0C(p*(e4+T“”' exp(-a?e4))-anQe^) (5.5.9) 

Bi = + GC(p*(e^/2+e^T"”*exp(-!!?e^))+ T^^expC-Te^) )-anQe^/2) 

(5.5.10) 

Substituting these values in equations (5.5,7) and (5.5.8), 
one gets 

dzn^ 1 'i 

= CG(p*(e+!D“ (exp(-Te) -exp(!De^))-anQ(e-e^)) 

(5.5.11) 

and o 

= QQ(pg("" Y" ^^'”^®xp(*“®«4) (e-e4)-^”'^ (exp(-Ie) -exp(-Ie^)'' 

“ (©•-64) ) + iSq (5.5.12) 
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5.5.3 P < e < 1 (Zone II) 


In this zone the circumferential stress is positive 
from P to 1, and the reinforcement provided in this zone 
is linearly vaiying from an^ to n^. So the equation 
(5.2.5) will be 

d^n^ 2 / = CG(pj(l-exp(-Te))-anQ-^hQ(e-p)) (5.5.13) 


First integration gives 

dm^2/‘^® “ 0C(p*(e+!D“^exp(-Ie))-anQe-EnQ(e^/2-pe)) + ^2 


(5.5.U) 

Second integration gives 

“x2 = CC(p*(eV2-T’'^exp(-Te)) 

+^ 2 © +B 2 (5-5.15) 

■vdiere 1-2 B 2 are constants of integration. 

From the continuity conditions (5.5.3) 

^2 = OC(p*(-e^-T”'‘exp(-Te^))-Hxnjje^-«;noP ^2) (5.5.16) 

B 2 = GO(p*(e^/2+!D"'’ e^ exp(*Te^) + !r'^exp( -!2e^)) 

^^ 4/2 + SiiQpVe) + m^j (5.5.17) 


Substituting these values in equations ( 5 . 5 . 14 ) stud 
(5.5.15) , one gets 



119 


dm^2 

de 


- CC(p*((e-e^) + (exp(-Te) ~exp(-Te^))) 


“Cxn^Ce— e^) — 11111^/2 (e—p ) } 


(5.5.18) 


and 


( ®~®4.) p 

m^2 ~ ^Pq ^ — 2 (exp(-Te) -exp(-Te^) 

-T""* exp(-!De^) (e-e^)) -anQ/2(e-e^)^ 

-Kno/6(e-P)^) + sIq (5.5.19) 

Using tlie boundary condition (5,5.5 ) in equation (5.5*19) 


= 


5ano(l-e^)^+Eno(l-P)^ 


== .,.11,1, I, , n . 1 „- n il, . I . M l,,, 

® (3(l_e^)2_6T-2(exp(-!r) -exp(-Te^))-6a?*’kp(-Te^)(l,.e^)) 

(5.5.20) 

5.5.4 e^iP (Case I) 

Por ^ P . Using fhe boundary and the continuity 
conditions (5.5.2) and (5.5.4) in (5.5.11) 


= 


°^^o^®5"*®4^ 


(6^-6^)'+ (exp(--Te 5 ) -exp(-3}e^) ) 


-1 


GG = 


2(mo+5o^ 


(anQ(e5-e^) ~p§((e5-e^)^-2T"**exp(-Te^) (e^-e^) 


-2a}’‘^(exp(-!De^) -exp (-Te^)))) 


(5.5.22) 
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when P ^ ^ 1 ■ 

Using the boundary and the continuity conditions (5.5.2) 
and ( 5 . 5 . 4 ) in equations (5.5*18) and (5.5.19), one gets 

- - ( 5 . 5 . 23 ) 




+ En^Ce^-P)' 


^((e^-e^) + T"’^(exp (-Te^) -exp (-Te^))) 


and 


GO = 


6(mjj+5Q) 


( 3 anQ ( 65 - 6 ^ ) ^+£no ( e^ -P ) ^ -p* ( 3 ( 65 - 6 ^)^ -60! e xp { -T p „ ) _ 


exp(-l?e^) ) "^^ 6xp(-!Pe^) (e^— s^))) 


( 5 . 5 . 24 ) 


5^5.5 P < (Case II) 


2 or this case, constants of integration A 2 and B 2 
for equations (5.5.14) and (5-. 5. 15) from boundai^y conditions 
(5.5.2) will be 

^2 ^ “C!G(p*(e^+T“"'exp(-Te^))-an^je^-KnQ(e^/2-Pe J) 

(5.5.28) 

and 

®2 " mo+CG(p*(e|/2+T“^exp(-Te^) + exp(-Te^) ) 

-^n^e|/2-Kn^j(e^/5-Pe|/2)) (5.5.29) 

Substituting these values in (5.5.14) and (5.5.J5) , one gets 


^2 


GG (p* (e+T”^ exp(-Te )-e^-T’'^ exp(-le^)) -an^ (e-e^) 
-Kng(eV2-Pe--e|/2+Pe^)) (5.5.30) 
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(e-e*)^ P i 

(Pq ^ — 2 ^ exp (-Te )-«xp (-Te 4) ) -!!!“' exp (-f (e-e^) ) 

-a2iy2(e-e^)^-Kno/6(e^-3PeJ3ee|-6pee^+2e|-3§e^) + 


( 5 . 5 . 31 ) 


Using tlie continuity conditions ( 5 . 5 . 5 ) 

3ano(l~e^)^+Knp(l-3P+6pe4.-3Pe^-3e4+ 2e| ) 


= 


( 3 ( 1 -e^) ( exp ( -T ) -exp ( -Te ^ ) ) -61"'’ exp ( -Te ^ ) ( 1 -e ^) ) 


(5.5.32) 

Using the continuity conditions ( 5 . 5 . 4 ) 

/2 + P(eg-e^)) 

((e^-e^) + (exp(-Te5)-exp(-Te^)) 




and 


CO = 


^(%+mo) 


(3anQ(e5-e^)^-KnQ(e|-3§e|-3e|+6Pe^-3Pe|+2e| ) 

-P^ ( 3 ( e^-e^) ^-6T”^ ( exp ( -Te^ ) - expt 5 e^ ) ) - 63 ?"^ exp ( -iTe ^) ( e^-e^} ) > 

( 5 . 5 . 54 ) 


5 . 5*6 Special Case 


let 


“o = “0 ■" “0 = ^ 


0 “O 

Therefore ^ P » equation (5.5.20) reduces to 


p-ji — 


3a(l-e^)^+ K (l-P)' 


(5(l-e^) -6T“^(exp(-T) -exp (-Te^))-6T”’ exp(-Te^) (l-e^) ) 

. ( 5 . 5 . 35 ) 

when e^ J P , equations (5.5.21) and (5.5.22) reduce to 
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p* = 
•^0 




(5.5.56) 


(e^-e^) + a?-"* (exp(-Te^) - exp(-Te^)) 


and 


GO = 


(a(e5^e^)‘^-p*((e5-e^)^ -2T”^exp(-Te^)(e5-e^)-2E 


,-2 


_2T"^(exp( 2e^)- exp(-Te^)))) 


(5.5.37) 


when P ^ , equations (5.5.23) and (5.5.24) will be 

2a(e^-e^) + K (e^-P)^ 

/ \ 1 

2(e^-e^) + 21*' (exp (-Te^)- exp(-Te^)) 

(5.5.38) 

and 


CG 


12 

(3a(e^-e^)^+K(e^-P )^-p*(3(e^-e^)^-6!D“^(exp(-Te5)-exp(!De^) ) 


-63'"'^exp(-Te^) (e^-e^))) (5.5.39) 

when P < e^ , equations (5.5.32), (5.5.35) and (5.5.34) 
will be 

3a(l-e^)^+K(l-3P+6Pe^-3e^ + 2e^ - 3Pe| ) 

3 ( 1 -e^) ^-61"^ ( exp ( -T ) -exp (-Te^)-6'D’”^ (exp( -Te^) ) ( 1 -e^) 

(5.5.40) 

2oc(e.-e.) + K (el - e^)-2Kp (e^-e.) 

p# = — 2— 3- ^ <5.5.41) 

(e^-e^) + 0!“”* (exp(-a!e 5 ) -exp (-Te^)) 


and 
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.2 


CC = 


(3aCe5-e^)^+K(e|-3Pe5-ro^+6N^-3Pe4+ 2e|) 
-V^(3ie^-e^)-ST~^(exv(-’Se^) -exp (-Te^)) 
-6T"”*exp(-Te^) (e^-e^))) 


(5.5.42) 


^en a = p = 1 


V * 


■■ A " " 

Cl-e^) -2T (exp(-T) -exp(-Te^)-2T“ exp(-!Oe^) (l-e^) 

(5.5.43) 


(1-64)' 






(e^-e^) + T“^ (exp (-Te^) - exp (-Te^)) 


(5.5.49) 


CG = 


((e^-e^)‘^-p*((e^-e^) -2T”^(exp(-Te^ ) -exp (-Te^)) 


-2T"^exp(-Te^) ( -^-e^))) 


(5.5.45) 


5.6 CONGLUSIOIT 

Three graphs have been drawn between CG and p* 
for different values of T, a and P. Pig. 5.6.1 is drawn 
for oc = P = 1 and for different T, It is clear from this 
graph that as the value of T in'- -cases ^the range of first 
mode increases. When the critical points of all the T 
are jot^aed for a = P = 1, one gets a ciarve #iich seperate 
the mode I from mode II and the equation of this curve ns 






p* Vs CC tor different T at a{. = p = 
for Janssen's pressure distribution 
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500,6416 p 

P* = - 55-3041 + 0.8929 OCR + - 0.0072 OCR 

° OCR 

( 5 - 6 . 1 ) 

In Pig. 5.6,2 for a = 0.75 and P =0,5 one can see 
that there is sudden drop in pressure in mode II at around 
00 = 76. The third mode for this starts at 288 and there 
is sudden rise in p*. So this curve is instable between 
00 = 76 to 288, Also one can see that the curve a = 0.5 and 
P= 0,25 and oc=p =0.5 cuts each other at CO equal to 510 and 
after that a = 0.5 and P = 0.5 gives more pressure. 

In Pig. 5 -6. 3, there is sudden drop of pressure 
in second mode of failure for a = 0,75 and p = 0,5 at 
00 = 48 and the third mode starts from 93 and there is 
sudden rise in pressure. So this curve is unstable 
in this range. 

This sudden drop in pressure is due to the 
reason that the value of ep lidiich is less than P before 
the drop in pressure becomes greater than p after the 
drop. As the value of T increases, the range of second 
mode for a = 0.75 and P = O .5 is shrinking. 
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Vs CC for variable 

*■^0 

for Janssen's pressure 







6. LIMIT LBSIGN OF SILO 


6*1 GENERIL 

A structure is to be designed for a satisfactory 
performance under different forces acting on it. It should 
not only serve the purpose for vdiich it is designed but 
it should also provide adequate safety against failure under 
most severe load conditions. The limit state design 
considers different ’failures’ at service conditions and 
collapse states besides dizrability^^^\ 

In this thesis only limit state design for 
strength has been considered. The working stress design 
is available readily in various standard bookS'^^’"*^^. 

The pressure distrioution in silos depends upon 
the type of material to be filled, type of silo walls, 
emptying condition etc. The pressure in silo during 
emptying is larger than the pressure during filling'' ^ . 

This increase in pressure depends upon the type of flow 
during emptying. There are generally two types of flow, 
non-d 3 nigmic and dynamic, the second flow gives larger 
pressure increase than the first flow. These flov® depend 
on Aether the material stored is compacted or loose prior 
to emptying* If the material before emptying is compacted. 
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it will give raise to nond 3 maiiiic flow and if it is loose, 
then the flow will be dynamic, There is no method to 
prevent dynamic flow. 

Here the Janssen's theory for pressure distribution 
in silos has been considered for design, ^ich gives 

(4) 

pressure at rest or for nondynamic type of flow 

Here the results from chapter five have been 
simplified for limit-state design. Only the value of 
Mp and Hp will change and so the value of GO, 

= K d^ f^j^ (6.1.1) 

= 0.87 -^g^fy = 0,87 M- 0 -t; t fy , for tension 

( 6 . 1 . 2 ) 

= 0*4- f^ij. t(l- Qj^) + 0.67 fy , for compression 

(6.1.5) 

T«h.ere 

f^^ = characteristic compressive strength of concrete, 

fy = characteristic yield s-trength of steel, 

d = effective thickness of shell wall, 

t == total thick. less of shell wall, 

^ Q-{j — ratio of steel to concrete for circumferential 

tension, 

P 0 Q = ratio of steel to concrete for circumferential 

CQBipression, 
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K = depends upon the type of steel and concrete. 


It has been assumed in analysis that = ^ 0 ^ = Ji 0 » 

to simplify the equation and also because of the reason 
that the compression zone is very anall and confined to 
0,1 Ii of the upper portion. 


Now considering the equilibrium equation 




= 1.5 P 


(6.1.4) 


where 1.5 is partial safety factor for static condition. 


Now 


M. 




^0 = 


M. 


rb 


9 




0 

pR 


N, 


(6.1.5) 


Substituting these values in equation ( 6 , 1 . 4 ) 

» 

■ck + nQ(0.87^iet fy)= 1 . 5 p*(O. 87 li 0 t fy) 


Kd^ f^v R 


de^ 


M iolc II 




R^(0.87 Pgt f_) ^e‘ 


+ Uq = 1.5P^ 
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— ^ + CO iiq = 1.5 CG p* (6.1.6) 

de^ 

where 

CO (°-87 l^fl t f^) 

Kd® f R 

ajad e = ■'■■ " ' (6.1*7) 

Jj 

6.2 NORIL4LISATION QP JAITSSEN'S EQUATIONS AND YIELD 
C01®ITI0N FOR SILO 


6.2.1 Bquilibriim Bq.uation 


Consider the equilibrium equation (6.1.6) 

d^ 

— r2-~. + GO ne = 1.5 CC p* (6.1.6) 

de2 ° 


where p 

2 Y 

p* = (l-exP(rTe)), 

PV (0.87 P gt fy ) 


T = P' k L /R^ 
and _ 

p’ = coefifieient of friction between stored material 
and silo wall, 

k = ratio of horizontal static pressure to vertical 
static pressure due to stored material, 
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= hydraulic radius of the silo 
Y = unit weight of stored material 
Row 

P* = 1.5 Pj (1 - exp(-Te)) (6.2.1) 

where 

P* = 1.5 P* = 3YRf/ (0.87Pet £y) 

So eq^uation (6.1.6) will be 

^2 
^ mx 

— + CC Uq = CO p* (l-exp(-Te)) (6.2.2) 

de'^ 

All the other equations and results of chapter 
five will remain the same for all modes because the initial 
equation is same. Here only the results for the case of 
<^ = P =1 have been listed for the design purpose. 

6.2,2 Results for Mode I 

Ydien m^ = m^ = = 1 and n^ = 

m^l = CG(p:»(eV2-T"^exp(-!re)-T"'‘e+T-2)-anQeV2) (6.2.3) 

p* = (6/OC + 5uQ)/(3(l+2I”^(l-exp(-T))-2T”'') ' (6.2.4) 

Pu = %®o (6.2.5) 

and 

OCR = 2/(nQeQ-p*(eQ+2T“^(l-exp(-TeQ))-2T"^eQ) (6.2.6) 
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6,2.3 Results for Mode II 


When = hIq = = 1 and = p iJ-^j 

For e ^ e-j 

nixi = GG(p*(eV2+T"^(l-exp(-Te))-T“‘'e)+nQeV2) 
For e 2 6-j 

m^2 = GG (p*(eV2+T"^(l-exp(-Te))-T"^e)-nQeV2 

.+(ee^-e^/2) (l+n^)) 

_ 2/QQ nQ-(l-HnQ)(2e.^-ef) 

" (l+2T"'^(l-exp(-T))-2T''‘ ) 


( 6 , 
(6 . 


n. 


K 


>62 - (1+^0^ 


62 + 5 !’"^ (exp(-Te)-l ) 


( 6 , 


CO = 


(Uoe2-(l+no)(2e^e2--ep-p*(e|-2T”'^exp(-Te2) 

(6. 




CCR=: 


-2!r“’*e2+2T“^)) 

iiQ(e2-e^) 


"1 

((e2--e^) + T” (exp(-Te2)-exp(-Te^) ) 

2 

(p*{e5+2T"^(l-exp(-Te3))-2F“'‘e5)-nQe5 
+(l+nQ) (26^63-6^)) 


( 6 . 


( 6 . 


(6.2.7) 

2 . 8 ) 

2.9) 

2 . 10 ) 

2.11) 

2 . 12 ) 

2.13) 
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6,2.4 Results for Mode III 


V/lien = 1 and = p M- 

(b e 

= CG(p*( — — -T“^exp(-0?e^)(e-e^)“l~^(expC-Te)-exp(-Te^)> 


p* = 


-nQ/2 (e-e^)^) + 1 

nn(l -e^)^ 


V* = 


-^C-'expC-Te^Xl-e^)) 
^o(^5 - ^4 ) 


^ ®5“® 4 ( exp ( -2?e5 ) -exp ( -T ) ) 


(6,2.14) 


( (l-e^)^-2T‘"^(exp(-T)-exp(-I!e^) ) 

,-1 


(6.2.15) 

(6.2.16) 


4 

CO = — — — 

(^o^®5~^4^ “P*( (e^~e^) — 2T exp(— !P6^)(e^“e^) 

-2T“^(exp(-Te5)-exp(-Te^))) (6.2.17) 

6,2,5 Yield Oondition 


The yield condition for the case of silos is 

m^(l+a)+2n^+2n^(2pp^+a-l)+2P^^^-4 pp^-2a = 0, for > 0 
and 

-m^(l+a)+2i3|+2n^(a_l)_2a=0 for < 0. 

This cylinder is bounded by the planes with equations. 

u-0 = 1 for compression 
hg == - d - PPq for tension , 
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where 


N, 




for GOTTipression 


°-«okt(l-4£o) + 0-6'7l" xc ^ 






% = 


= 


^6 = 


o^erju^^ f 


for tension 


M. 


Kd f 


ck 


e 


^y 

N, 


for tension. 


0.+f„i.t(1-MB„)+0.67tiBot 

Assume that 


for compression , 


= Ajj/t , reinf orcement ratio in axial direction, 

^ec = ^et = = Ag/t , reinforcement ratio in circumferential 

direction, 

a - tensile yield strength of concrete 

compressive yield strength of concrete 

yield strength of reinforcement 

p = ■ ■ - " • 

compressive yield strength of concrete 

6,2,6 Interaction Equations 

The interaction equations for combined bending 
and axial force for balance force pair, i.e. the section . 
f ails simultaneously in compression and in tension, are 
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^ab 

= 0.36 S^bd + 0.87 Usc-^3t^ 

(6.2.18) 


^ab 

V Kb +0.87 A3^ fy 

(6.2.19) 


^ao 

= (1-Px) (0-45 b t f^j^) + 0.75 fy btp^ 

(6.2.20) 


®x 

II 



P^ 

^ao 

(6.2.21) 


3p 



If 


C i 4- ^ 

< Pg_ then the design is adequate^ 


where 




= 

distance of Pg^^ from tension steel , 



= 

0.5 Bg + ejj , 


^s 

= 

distance between compression and tension 

steel, 

®b 

= 

distance of P^^ from centre of the shell 

wall, 

■^SC 

= 

area of steel in compression, 


-^st 

- 

area of steel in tension, 



= 

depends upon the t 3 rpe of steel, 


®x 

= 

actual eccentricity, 



= 

ratio of area of asiel reinforcement to concrete, 


= actual axial capacity of the silo wall along with 
moment , 

Pg^Q = ajtial load capacity with zero eccentricity, 

Pg^j^ = axial load capacity with eccentricity 
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6,3 DESIGN. OB' SILO 
6.3.1 General 

Here, design is based on static pressure distribution 
only. No correction factor lias been considered for wind 


force, e art hqualce force 

or emptying condition. 

DATA: 

Material stored 


Cement 

Concrete mix 

• 

M20 

Steel 

• 

* 

Grade Re 415 

Unit weight of cement 


U.l IcN/m^ 

Ratio of horizontal static 


pressure to vertical static 
Pressure due to stored cement : 

0.27 

Coefficient of friction 

betv/een 


cement and silo wall 

« 

0.70 

Minimum reinforcement 


0,25 percent 

Short silo Diameter 

: 1.5m 

Height ; 5 m 

Medium silo Diameter 

: 3,0m 

Height :10 m 

Deep silo Diameter 

: 6.0m 

Height :2C m 


Axial pressure due to friction is 
e 

= D /YR^d^expC-Ie)) de = YR-,L(e+!D“^ (exp{-Te)-t )) 

(6,3.1) 


#iere 1 = total depth of silo 
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So 

T = 2.52 P 

I^gt 

CC = 130.815 ( — — p -) 

a = 0.0 
P = 20,75 

Now the yield condition will be 

2n.^+2n^iA-^,5\^^-^)+ 861 . 125 m ^ -83\x^+m^ = 0 

for lOjj > 0 ( 6 , 5 . 2 ) 

2^x"2^"“x = 0 f or < 0 (6.3.3) 

and ng = 1 fOr compression 

ng = - 20.75 Hq for tension 

Moment distribution for first mode, second mode 
and third mode have been shown in Fig. 6.5.1 for T = 2.52 
and m^ = 5^ = n^ = Eq = 1 for GC = 20, I 30 and 525 respectively, 

6 , 3,2 Design of Short Silo 
6 , 5.2,1 Design of wall 
Let , 

Percentage of circumferential reinforcement : 0.25 

Total thickness of silo wall 0,12 m 

Effective t;^ckness of silo wall i 0.10 m 

So 

CC = 150.815 

and T = 2.52 








M 




HO 

By using equations (6.2.3) through. (6.2.6) 

P* = 0.1354 

and p* = 3Y R^/ (0.87 b fy ji') 

Prom this P’0 = 0.1 45 percent. 

Minimum required rei’iforcement is 0,25 percent and 
provided is also the same. So circumferential reinforcement 
is adequate. Now for all m^ , is calculated for known 
n0 from equations (6.3.2) and (6,3.3), fhe required n^ 
is calculated from equations (6.1.3) and (6.3 .I). All 
these stress resultants n0, and also required 

for this case have been presented in Pig. 6.3.2 for 

“ '•.75 percent and p0 = 0,25 percent in the complete 
silo. 

Now, = 1.5 K b d^ ^ck “ ^'••4- kNm/m 
where , K =0.138 for the given steel. 

Por this case, critical situation is at bottom 
of the silo because of the maximum positive moment and 
the maximum axial force. 

At e = 1 

Nx = 47.7 kN/m 
and = 44.4 ld!Tm/m 

is positive at this level so n^ calculated from yield 
condition as well as n^j required will depend upon P^j. 
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Provided = 1.75 percent at bottom, out of 
which = 1,25 percent and = 0.5 percent. So from 
equations (6,2.18) through (6.2.21). 

= 88.911 k!T/m 

^ab 4-4.93 kNm/m 

= 0.465 m 

e_ = 0,868 m 

Pg^ = 48.88 Idl/m > 47.7 IcN’/m so safe. 

The other critical point is at e =0.66, Here 
the calculated is less than the required n^.. The 
moment is changing from positive to negative. If the 
is changed , calciilated will not change except n^ 
required. Here, 

= 0.013175 = 0.545 kNm/m 

= 28.27 kJf/m 

Alternate bars from outer face and alternate two bars 
from inner face have been curtailed from e = 0.75 onward. So, 
^xt “ *^*^5 percent and = O.4I66 percent. 

How from equations (6.2.18) throu^ (6.2.21) 

^ab = 486.037 Idf/m 
^ab®s= 35.139 Mm/m 

e^ = 0.0325 m 

^ao = 1356.42 M/m 
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e„ = 0,0193 m 

Pg^ = 655.44ld?r/m > 28,27 Idl/m so safe, 

Now from e = 0.5 onward alternate bars from 
inner face have been curtailed, so the reinforcement at 
e = 0,43 where the maximum negative moment occurs is 
M' = 0,25 percent and ~ 0,2083 percent, 

= 0.167 = 6.92 kNm/m and = 16.315 Wm. 

Prom equations (6,2. is) through (6.2.21), 

P^ = 26.7 kN/m > 16.315 kN/m. So safe 

3L 

Reinforcement details: 

Provide 8 0 bars at 300 C/C as circumferential 
reinforcement on both faces. Provide 10 0 bars at 50 G/G 
as axial reinforcement on inner face and at 130 C/G outer 
face from e = 1,0 to 0,75, Provide 10 ^ bars at 150 C/G 
on inner face and 260 C/C on outer face from e = 0.'^5 to 0.5. 
Prom e = 0,5 onward provide 10 0 bars at 300 C/C on inner 
face and 260 C/C on outer face. 

6,3* 2,2 Design of conical dome and ring beam 
Data: 

Diameter of opening provided in conical dome : 0.25 m 

Angle of inclination of conical dome wall 

from vertical axis : 45° 


fhickness of the wall 


0.12 m 
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= 47.7 Wm 

= P^ tan e = 47.7 kN/m 
H 2 = 1.5p = 10,425 

vfcere 

P^ = vertical force due to wall and stored cement, 

= radial force due to the action of vertical force> 

H 2 = radial cement pressure, 

0 = angle of inclination of cone from vertical axis, and 

b = width of ring beam. 

Let b = 0,2 m 

Total hoop tension = D/2 (H^ + H 2 b) = 37 kPT 

2 

So steel required in ring beam = 102,5 mm 

Provide two 10 0 bars in ring beam excluding circumferential 

reinforcement. 

Slanted length of the conical dome = 0,884 m 

Diameter at mid hei^t = 0,875 m 

Wg(weight of conical hopper) = 7.3 MT 

WjjjCwei^t of material stored in conical hopper) = 6,2 kfT 

^des = — F“ ^ ^ 

where ~ vertical static pressure at the top of the 

hopper. Substituting the values, one gets 

qaes = 25.75 H/m^ 
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P 


mu 


= 1.5 ( 




D 


4 sin 8 


W, 


m 


W, 


7iD sin 6 


tcD sin 9 


) 


where, = diagonal force per unit length, and D = diameter 
of the hopper at the level being considered. 


^mu ^ hopper ) = 26.56 kN 
^mu middle of the hopper) = 22.32 kN 
^mu. bottom of the hopper) = 39.88 kN. 


So the required meridioxaL reinforcement at bottom, 

2 

A„ = 1 10.5 mm /m, 
s 

This is less than the minimum reinforcement, so 
provide minimum reinforcement i.e, 0,25 percent. So 
provide 8 0 bars at 165 G/C at outer face only. Curtail 
alternate bars from the point vdiere the diameter reduces to 
0.75 m. 


^tu 


qde3° 
2 sin 0 


where = hoop force per unit width. 

Substituting the values, one gets 


P-fcu =41 ijN/m 

2 

So the required hoop reinforcement, Ag = 113.47 mm /m^ 

This is less than the required minimua, so provide 
minimum required reinforcaaent i.e. 0.25 percent. Provide 
8 0 bars at 165 C/G at outer face only. 

The reinf credent details are showi in Pig, 6,3.5. 
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6, 5- 2, 3 Cost computation 


Unit costs of the materials including fixing, or 

(14) 


placing or finishing etc. are 


One m of M20 concrete 


= Rs. 450.00 


One m of steel 
Volume of concrete in wall and ring beam = 2.928 m 


= Rs. 3554.0 
3 


=0.292 m' 


= 3.22 m-' 

= Rs. 1450.00 
= 0.0321 m^ 

= 0.00161 
= 0.03371 
= Rs. 1131.00 
= Rs. 2581.00 
= SI .-27 5 


Volume of concrete in hopper 
Total volume of concrete 
Cost of concrete 

Volume of steel on wall and ring beam 
Volume of steel in hopper 
Total volume of steel 
Cost of steel 
Total cost 
Total stored cement 
Cost of storing for one cubic meter of cement = Rs. 278,28. 

6.3.3 Design of Medium Silo 

6,3*3. 1 Design of wall 

let , 

Percentage of circiaaferential reinf oreeanent : 0,5 

Total thickness of silo wall : 0.12 m 

Effective thickness of silo wall : 0,10 ® 

So, CG = 523.26 and T = 2.52. 
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By using equations ( 6 . 2 , 7 ) through (6.2,13) 

= 0.18455, so>Q = 0.43 percent. Provided is 

0.5 percent, so the circxomf erential reinforcement is 
adequate. Now for all m^, n^ is calculated for known ng 
from equations (6.3.2) and (6.3.5). Also the required n^ 
is calculated from equations (6.1.5) and (6.3-1). All these 
stress resultants n0, m^, n^ and also required n^ for this 
case have been shown in Pig. 6,3.5 for ^0=0. 5 percent 
and ~ percent throughout. 

= 1.5 K b ^m/m. 

Por this case, critical situation is at bottom 
of the silo because of the maximum positive moment and 
the maximum axial force. Now at e = ^ , N^ = ^ 4 . 5 ,1 
= 41 ,4 kNm/m. 

The reinforcement provided on inner face is 1.25 percent 

and on outer face is 0.5 percent. So ^^xt “ 1*25 percent 

g tr id P = 0.5 percent . So from equations (6.2.18) through 
xc 

(6.2.21). 

Pq == 146 kN/m > 145.7 kN/m. So safe. 

A.t e = 0.89, n^ calculated from yield condition is smaller 
than the required. This is because moment is ohanging from 
positive to negative. Here the tension face is becoming 
compression face and compression face is becoming tension 
face so , =3 0,5 percent, P^jQ ^ 1.25 percent. 










wm§. 




/'V/4^ 






Now from equations (6,2,18) through (6,2,21 ) 

= 1628 li2l/m > 126,8 kN/m , so safe. 

At e = 0.61 the negative moment is maximum i.e. 
and is equal to 74.8 ktT/m. Prom equations (6.2.18) 
through (6,2.21) 

P = 80.6 IdT/m > 74.8 kN/m. 

gL 

Here percent and = 0.625 percent. 

This is because from e = 0.8 alternate bars from inner face 
have been curtailed and alternate bars have been provided 
on outer fane. 

At e = 0,21 the n^ calculated is less than the n^ 
required. Here ^xc ~ 8.3125 percent and ^xt = 0.5 percent 
because alternate bars from both face have been curtailed 
from e = 0,5 onward. Prom equations (6,2.18) throu^ (6,2.2l) 

P =144.3 ktf/m > 17 PN/m. So safe. 

51 

How curtail alternate bars from outer face from 
e = 0,15 onward. At e = 0,1 there is relative maximum 
moment. Prom equations (6.2.18) through (6.2.21), one gets 

!*_ = 32.9 PH/m > 6.4 PH/m, So safe 

d 

Reinforcement details: ' 

Provide 8 0 bars 165 C/C as circumferential reinforcement 
on both faces. Provide 10 bars at 50 G/G as axial 
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reinforcement on inner face and at I 3 O C/G on outer face 
from e = 1.0 to 0,8, Provide 10 bars at 100 C/C on 
inner face and 65 C/C on outer face from e = 0.8 to 0.3, 
Provide 10 0 bars at 200 G/C on inner face and I 3 O C/C 
on outer face from e = 0.5 to 0,15. Prom e = 0,15 to top 
provide 10 0 bars at 200 C/C on inner face and 260 C/G on 
outer face. 

6, 3 . 3,2 Design of conical dome and ring bean 
Data: 

Diameter of opening provided in conical dome : 0. 35 m 

Angle of inclination of conical dome wall from 
vertical axis : 45° 

Ihickness of the wall : 0,12 m 

P^ = 145.75 M/m 

= P^ tan e = 145.75 M/m 

H 2 = 1.5 P = 20,8 M/m^ 

Let b = 0,2 m ' 

Total hoop tension = 225 M 

2 

So steel required in ring beam = 623,2 mm 

Provide six 12 0 bars in ring beam excluding circumferential 
reinforcement. 

Slanted length of the conical dome = 1.88 m 
Diameter at mid hex ^t = 1.675 m 

Wg (wei^t of conical hopper ) =29,7 kH 

Wjjj (wei^t of material stored in conical hopper) = 50 M. 



151 


laes = 51.5 tN/m^ 

hopper) =100 kH/m 

^mu middle of the hopper)= 79.9 kN/m 

^mu bottom of the hopper )= 165.26 M/m 

Reinforcement required at top of the hopper = 276,97 mm /m. 

This is less thatn the minimum requirement, so 
provide minimxim requirement of steel i.e. 10 JZi bars at 
260 C/C on outer face only. Reinforcement required at 
bottom of the hopper is 452,2 mm /m. So needed reinforcement 
is 10 0 bars at 170 C/C on outer face only. Now provide 
10 0 bars at 260 0/C at top of the hopper. Cuirfcail alternate 
bars from the place vfeere diameter of hopper .reduces to 2 m 
and again alternate bars from the place vdiere diameter 
becomes 1 m. So spacing of 10 0 bars at bottom is 120 0/ 0. 

Now = 164 kN/m , so the required hoop reinforcement 
is 453.95 mm^/m. Provide 10 0 bars at 17 0 O/O at outer face 
as hoop reinforcement. For reinforcement details see 
Fig. 6.3.6. 

6. 3.3* 5 Cost computation 

Yoltane of concrete in wall and ring beam = 11,602 

■ ' 3 ' 

Volume of concrete in hopper = 1,188 m 

Total volume of concrete =12.79 m^ 

Cost of concrete = Rs.5756 
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Volume of steel in wall and ring beam 

Volume of steel in hopper 

Total volume of steel 

Cost of steel 

Total cost 

Total stored cement 

Cost of storing for one cubic meter 
of cement 


0.2232 m^ 
0.009504 m^ 
0,2327 

Rs. 7805 
Rs. 10561 
74.232 m^ 

Rs. 142.27 


6 , 5.4 Design of A Deep Silo 
6.3.4 . 1 Design of wall 
Let, 

Percentage of circumferential reinforcement : 1.0 

Total thickness of silo wall : 0.12 m 

Effective thickness of silo wall : 0,10 m 

So GO = 2093. O4 . By using equations ( 6 , 2 . 14 ) through 
(6,2,17), one gets 

p* = 0.3205, so P'0 = 0,979 percent. Provided is 

1 percent, so the circumferential reinforcement provided is 
adequate. Now for all m^ n^ is calculated for known n0 
from equations (6.3.2) and ( 6 , 3 . 3 ). The reqiiired; n^. Is 
calculated from equations ( 6 . 1 , 3 ) and (6 . 3 . 1 ). All these 
stress resultants n^* n^ and also required n^ for this 
case have been shown in Pig. 6,3,4 for fts = 2 perc^t and 
^0 = 1 percent in the complete silo. 
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At e = 1 , Hjj. = 496.6 kN/m and = 41.4 Idlm/m, 

Provided ^ percent at bottom out of vSiicb. Pxc ~ f^xt 

= 1 percent. So from equations (6.2.18) throui^ (6,2.21). 

Pa = 508.7 IdVm > 496.6 kN/m. So safe. 

Curtail one bar after every alternate three bars 
from both faces from e = 0.8 onward. At e = 0.75 moment 
is majcimum negative i.e. So from equation (6, 2. '18) 

through (6.2.21) for = 0.75 percent, 

P^ = 425.5 kS/m > 332 kU/m. So safe. 

21 

At e = 0.3 moment is maximum positive i.e. 
the same reinforcement. So from equations (6,2.18) to 
( 6 , 2 . 21 ). 

Pa = 134.2 kif/M > 83.36 k^T/m. So safe. 

Continue this same reinforcement upto top. 
Reinforcement details: 

Provide 10 0 bars, I30 C/C as circumferential 
reinforcement on both faces. Provide 10 0 bars 65 C/C 
as axial reinforcement on both faces from e= 1,0 to 0.8, 
Provide 10 0 bars 85 C/C on both faces from e = 0.8 to top, 

6,3,4,2 ^Design of conical dome and ring beam 

Riveter of opening provided in conical dotee 4 


Bata: 
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Angle of inclination of conical dome wall 

o 

from vertical axis : 45 

Thickness of wall : 0,12 m 

= 496.6 Idl/m 

= 496.6 IdJT/m 

Hg = 41.7 

let , b = 0.30 m. 

Total hoop tension = 1516 kN' 

2 

So steel required in ring beam = 4200 ram /m. 

Provide 9 bars of 25 0 in ring beam excluding 
circumferential reinforcement. 

Slanted length of the conical dome = 3.89 m 


Diameter at mid height 

= 3.25 m 

Wg - 

119.2 kN 


\ = 

398.5 kN 


9des~ 

103 kff 


Pjjj^Cat top of the hopper) 

= 386 IrN 


middle of the hopper) 

= 285 kN 


bottom of the hopper) 

= 726.5 kN 


2 '■ 

Reinforcement required at top is 1069,1 ma/m. So 
provide 12 0 bars at 100 C/C. Reinforcement required at 
bottom is 2012,2 

Provide 12 0 bars at 100 C/C at top. Curtail alternate 
bars f re® the place Miere diameter becomes ^ and again 
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cuartail alternate bars from the point vfcere diameter reduces to 
2 m. At bottom spacing is 34 C/C 4 So safe. 

= 655.5 MT/m , so hoop reinforcement required is 
2 

1815,5 mm /m. So provide 12 0 bars at 120 C/C at outer 
face as hoop reinforcement. For reinforcement details 
see Fig. 6.3.7. 

6.3.4 . 3 Cost computation 


Volume of concrete in wall and ring beam 

= 46.398 m^ 

Volume of concrete in hopper 

= 4.768 m^ 

Total volume of concrete 

= 51.166 m^ 

Cost of concrete 

= Rs. 23025 

Volume of steel in wall and ring beam 

= 1.28962 m^ 

VolumS of steel in hopper 

= 0.07947 m^ 

Total volume of steel 

= 1 . 36909 m^ 

Cost of steel 

= Rs. 45920 

Total cost 

= Rs. 68945 

Total stored cement 

=593.7484 m' 


Cost of storing for one cubic meter of cement =ae, II 6.124 




I 50c/c i:n',bot^!q<ti':200€ 
100 c/c irt nex* 5C00 
200c/c tn top 3000 - 

inrifir* ^ bofs , , , 

1U00U ^30c/c in bottom 200< 

I floSc'ln^pext fsO^O 

■■"I z&0c/c,»n iopiapo 

t , O'bbrs'oi;',:^-. 
200 . 




Im) I 


8 ^ at T 
165 c/c i 

\ i2ot: 






'fi^..6*3'5'"'-R'ew'f or cement details 


'0;S ;C;''/b''in bottom 4000 
iri"top;16000 ■ . 




12 

.120 c/c 
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6.4 DISCUSSION 

As the size of silo increases, the mode of failure 
changes from first mode to second mode or second mode to 
third mode. Also the cost of storing material per unit 
volume decreases with increase in the size of the silo. 
This is because of the better utilization of the material 
stren^h over a larger portion of the structure. 



7. DISCUSSION OB' RESULTS AND CONCLUSION 


7.1 GENERAL 

The basic differences between a silo and a tank 
are due to the nature of normal pressure distribution 
and the wall friction forces. In the case of a tank, 
storing fluid, the lateral pressure increases linearly 
with depth, v\iiereas it is nonlinear in the case of a silo. In 
the case of the former the rate of increase in pressure 
with depth is constant vhereas in the latter the rate of 
increase in pressure increases with depth. In the case 
of a fluid storage tank axial forces on the wall will be 
present only v\hen the fluid is hi^ly viscous, whereas in 
the case of silo axial forces due to friction are present 
for all materials like cement, coal, grain etc. In the 
case of fluid storage tank the difference between loading 
(filling) and unloading (exptying) pressure is negligible 
while in the case of silos the pressure on the wall during 
emptying process may be as hi^ as four times that of the 
fillingin pressure. This pressure increase during emptying 
will depend ant only upon material but also the conditions, 
namely, compact or loose packing^. The silo walls as siaaed 
to be free at the top and clsaaped at the bottom. However, 
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in practice, this is not 't;rue. 

7.2 DISGUSSOT OP RESULTS 

In all the storage structures considered herein, 
there are three basic modes of failure. A typical 
circumferential moment capacity distribution for a 
cylindrical wall will consist of three zones. Zone I 
is upto a depth e-jL therein the nature of hoop force is 
compression in Zones II and III there is hoop tension. In Zones 
I and II the hoop force is uniform vdiile in Zone III 
it is assumed to be varying linearly. The region boundary 
between the Zones II and III is pL from the top of 
the cylindrical shell. The typical modes of failure are 
showti in Eig. 3.3.1, wherein the dots indicate the 
location of the plastic hinge circle. The number of 
hinge circles in failijre modes indicate the mode of 
failure. In modes one and three the hoop compression 
zone is absent. In the case of the third mode, the 
top most portion of the wall beyond the failure zone 
upto a depth of e^L, any stress distribution that is 

statically admissible may be assumed (Fig. 3.7.1). 

From equations (4,2.6), (4.5.2) , (5.1.1) and 
(5.2,2) one can readily see that Janssen's theory is giv*- 
iiSg a lower pressure than the Reimbert 's theory for 
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lower value of nondiraensional friction parameter (T). As 
the value of T increases the difference decreases and 
for T > . 2.6 the Janssen’s theory will give more 

pressure than Reimbert ’ s theory. This difference increases 
for T around 6.5 and after this , the difference starts 
decreasing and for T more than 20 this difference is 
almost negligible. The Reimbert ’s theory is found to give 
better correlation between theoretical and experimental 
data. However in this thesis only Janssen’s theory, that 
is widely adopted, has been considered for working out a 
design example for T = 2,52 at Yhich the difference in 
pressure between these tv/o theories is almost negligible. 

So for smaller value of T the Reimbert ’s theory -is 
conservative when compared with the Janssen’s theory but 
for large value of T this is not the case, 

Rrom Fig, 3.8,1 it is clear that 62 e^ are 
always greater than P while e^ and e^ are always less 
than P, This is tmie for P _< 0.5 and a upto 0.75, Hence, 
these curves are smooth and changes in slopes are noticed 
only at critical points. It is also clear that the hoop 
compression is confined only in the top 10 percent depth 
of the tank wall. 

From Figs. 4 . 7.2 and 4.7,3 it is clear that e'|»®2>®4 
and e^ are always less than p for P < 0.5 and a 0,75. 
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Hence, these curves are smooth and changes in slopes are 
noticed only at critical points. In Fig. 4.7.2 the curves 
fora=p=I =1 and a = 0.75 , P = 0.5 and T = 1 intersect 
each other at the nondimensional shell parameter GO = 280. 
This is due to the different failure modes. For the first 
curve the failure mode is the third one while for the latter 
it is in second mode. 

In Figs. 5 . 6.2 and 5.6,3 it has been found that 

for the case of a = 0,75 and P =0.5 for T = 1 and 2.5 

for second mode of failure there is a sudden drop 'in 

nondimensional pressure parameter (?*)• I’rior to the drop 

the value of 62 is less than that of P and it is greater 

than p after the drop. As the value of nondimensional 

friction parameter (T) increases , the range of nondimensional 

shell pargpeter (OC) for second mode of failure decreases. 

Besides, vdien it approaches the third mode there is a sudden 

rise in nondimensional pressure parameter (P*). For T = 1 

at GO =76 there is a sudden drop in p* and for a value of 

o 

GO = 288 there is not only a jump in p* but also a chan^'B in 

o ^ 

the mode of failure. The corresponding values for 
T = 2,5 are OG = 48 and 93 respectively. 

In both the cases of silo it has been noticed that 
the values of e^, e^ , Sgj e^ and e^ decrease as this value 



163 


of T increases. This decrease in e 2 for the case of 
a = 0,5 and p = 0,25 is faster than the decrease in 
the rate of other cases, 

From section 6,3 it is clear that for a given L/R 
ratio the cost of storing per unit volume of material 
decreases as the silo beccmies deeper for oc = p = 1 , For 
L/D = 3.33, the unit cost for the three types of silos 
(shallow, medium and deep) are Rs. 278.28, Rs. 142.27 and 
Rs. 116,12 respectively. 

On comparing the results shown in Tables 4.1 to 4.3 
and 5.1 to 5.3 it is obvious that the critical values of CG 
(junction modes) for Reimbert’s pressure distribution 
theory are greater than that of Janssen's. As the value 
of T decreases the difference between the CGR values for 
the two theories will decrease. This can be inferred 
the forementioned tables. 

7 , 3 COIfCLUSION 

For a very deep silo both the Reimbert 's and 
Janssen's theories will yield almost the same results. 

For either of the pressure distribution,f or a given 

nondimensional shell parameter (GC) and nondimensional 
friction parameter (T), one can readily determine the 
mode of collapse by using the Figs. 4 . 7.1 to 4 . 7.3 and 
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5.6.1 to 5 . 6 . 5 . For a given size of tank one can 
readily determine the mode of collapse from Fig, 3 . 8 . 1 . 

The same set of figures mentioned above for 
predicting the failizre mode can be used to ' design a 
storage structure in conjunction with the yield 
condition described by equations(2.4.1 ) to (2.4,4), 
outlined in Section 6.5. 

The failure modes in the case of water tank 
and the silo are similar. For the first time, plastic 
analysis of silos is presented with variable circumferential 
reinforcement along the depth. Other conditions permiting, 
the deeper the silo the better it is in terms of economy 
(cost per unit volume of material stored), 

A shell of given dimension may fail in different 
modes for different pressure distribution theories. 

Hence, one should be ext3?emely careful in selecting a 
particular theory. Such a theory should mathematically 
represent the physical situation. 

7.4 SCOPE FOR FURTHER STUDY 

Silos covered at the top may be analysed. 

limit state design of the silo is the next ^ ^ 
logical step, wherein the collapse and the servioebility 



requirements are considered 


!Phe degree of fixity at the junction of the 
cylindrical silo wall and the conicsil hopper should 
be studied in depth. 
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